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Preface 


The  eutrophication/general  water  quality  aspect  of  die  New  York  Bight 
(NYB)  Study  was  conducted  at  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  for  the  U.S.  Army  Engineer  District,  New  York 
(CENAN).  The  study  was  funded  by  CENAN  under  Section  728  of  the  Water 
Resources  Act  of  1986.  Mr.  H.  Lee  Butter,  Chief,  Research  Division,  Coastal 
Engineering  Research  Center,  was  the  WES  study  manager.  Ms.  Lynn  M. 
Bocamazo  and  Mr.  Bryce  W.  WisemOter  of  CENAN  were  responsible  for 
overall  project  management  of  the  NYB  Study. 

This  study  was  conducted  by  Mr.  Ross  W.  Hall  of  the  Water  Quality  aid 
Contaminant  Modeling  Brandi  (WQCMB),  Environmental  Processes  and 
Effects  Division  (EPED),  Environmental  Laboratory  (EL),  WES.  Dr.  Mark  S. 
Dortch,  Chief,  WQCMB,  assisted  in  die  analysis  and  interpretation  of  results 
and  provided  direct  supervision.  General  supervision  was  provided  by 
Mr.  Donald  L.  Robey,  Chief,  EPED,  and  Dr.  lohn  Harrison,  Chief,  EL.  This 
report  was  written  by  Mr.  Hall  and  Dr.  Dortch.  Report  review  was  provided 
by  Mr.  Thomas  M.  Cote  and  Dr.  Carl  F.  Cerco  of  the  WQCMB. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robot  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


1  Introduction 


Background 

The  New  York  Bight  (NYB)  Study,  funded  through  Section  728  of  die 
Water  Resources  Act  of  1986,  is  an  investigation  of  die  technical  feasibility  of 
conducting  hydro-environmental  modeling  and  monitoring  of  the  NYB  system. 
The  modeling  efforts  of  this  study  were  focused  on  several  aspects,  including 
hydrodynamics  and  circulation,  eutrophication  and  general  witter  quality,  con¬ 
taminants  and  toxic  substances,  and  particle  tracking  for  dissolved  aid  sus¬ 
pended  matter,  such  as  floatables.  This  report  presents  the  results  of  the 
eutrophication/general  water  quality  model  study. 

In  general,  eutrophication  refers  to  increased  productivity  and  degraded 
health  of  a  water  body  as  a  result  of  excessive  nutrient  inputs.  The  health  of  a 
system  can  be  characterized  in  terms  of  the  diversity  and  abundance  of  biolog¬ 
ical  communities.  The  term  general  water  quality  is  used  here  to  refer  to  the 
more  conventional  water  quality  constituents  that  are  used  to  measure  the 
health  of  a  water  body,  such  as  dissolved  oxygen  (DO),  nutrients,  and  algae. 
The  NYB  is  not  considered  globally  eutrophic,  but  there  are  regions  that 
experience  low  dissolved  oxygen  in  near-bottom  waters,  which  is  detrimental 
to  marine  life.  There  is  legitimate  concern  that  human  activities  such  as 
municipal  and  industrial  wastewater,  combined  sewer  overflows,  stormwater, 
and  various  dumping  operations  may  be  increasing  die  eutrophication  of  the 
Bight,  as  has  happened  in  other  coastal/estuarine  systems  around  the  United 
States.  The  intent  of  die  Section  728  legislation  is  to  initiate  activities  that 
could  lead  to  improved  awareness  and  protection  of  the  New  York  Bight’s 
health.  Numerical  simulation  models  offer  die  most  cost-effective  and  techni¬ 
cally  defensible  means  of  assessing  die  impacts  of  natural  events  and  human 
activities. 

Section  728  mandated  studies  to  investigate  the  technical  feasibility  of 
developing  methods  (e.g.,  modeling  and  monitoring)  to  assess  environmental 
conditions  in  die  Bight  and  how  the  Bight  environment  will  respond  to  human 
influence.  The  term  “technical  feasibility”  is  used  because  the  NYB  system  is 
very  large  and  complex.  The  system  consists  of  tidally  influenced  estuaries, 
harbors,  and  bays;  Long  Island  Sound  (LIS);  the  Apex  region  between  the 
open  waters  of  the  Bight  and  die  harbore/estuaries;  and  the  Bight,  which 
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extends  (for  this  study)  from  Cape  May,  NJ,  northeasterly  approximately 
550  km  along  the  coastline  to  Nantucket  Island,  MA,  and  extends  approxi¬ 
mately  160  km  offshore  to  beyond  the  continental  shelf.  The  depth  of  the 
study  site  varies  from  3  m  to  over  2,000  m  seaward  of  the  continental  shelf. 
The  system  is  influenced  by  many  physical  processes,  such  as  astronomical 
tidal  fluctuations,  meteorological  forcings,  and  large-scale  oceanic  circulation 
patterns  affecting  the  entire  Middle  Atlantic  Bight  Significant  local  effects 
include  riverine  inflows  and  bathymetric  variations  (Scheflher  et  al.,  in  prepa¬ 
ration).  Water  quality  is  affected  by  the  above-mentioned  physical  processes 
plus  turbulent  diffusion,  material  exchange  with  the  atmosphere  and  bottom, 
and  chemical  and  biological  processes.  All  of  these  factors  and  processes 
present  a  substantial  challenge  for  simulating  the  NYB  system. 


Objective 

The  objective  of  this  study  was  to  determine  the  technical  feasibility  of 
developing  a  numerical  eutrophication/general  water  quality  simulation  model 
to  assess  the  impacts  of  natural  and  human  activities  on  the  NYB.  The  best 
way  to  investigate  the  feasibility  of  something  is  to  create  and  develop  it  and 
test  it  This  was  the  approach  taken  here. 


Approach 

The  modeling  technology  recently  developed  for  Chesapeake  Bay  (Cerco 
and  Cole  1989;  1992;  in  preparation)  was  applied  to  the  Bight.  This  technol¬ 
ogy  consisted  of  three-dimensional  (3-D),  time-varying  hydrodynamic  and 
water  quality  models.  The  hydrodynamic  model  (HM),  which  is  described  in  a 
separate  report  for  the  study  (Sche finer  et  al.,  in  preparation),  provides  the 
circulation  required  by  the  transport  terms  of  the  water  quality  model  (WQM). 
The  WQM  is  indirectly  linked  to  the  HM.  Thus,  the  HM  is  applied,  and  the 
output  is  stored  and  subsequently  used  as  input  for  the  WQM. 

A  workshop  on  Bight  modeling  was  held  in  New  York  City  (WMI  1989)  at 
the  beginning  of  the  Section  728  studies.  One  of  the  primary  recommenda¬ 
tions  of  this  workshop  was  that  the  models  should  include  the  harbors  and 
estuaries.  Long  Island  Sound,  and  the  Bight  out  to  the  shelf-break.  It  is 
important  to  capture  the  circulation  among  these  three  regions  (e.g..  Bight, 
Sound,  and  harbors/estuaries).  Therefore,  this  model  does  include  all  three 
components. 

Both  the  HM  and  WQM  are  mechanistic,  deterministic,  numerical  models 
based  on  conservation  principles.  The  WQM  is  based  on  foe  conservation  of 
mass  and  includes  sources  and  sinks  of  mass  arising  horn  kinetic  reactions, 
transfers,  and  transformations. 
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Application  of  the  modeling  technology  to  the  Bight  proceeded  through  the 
following  steps: 

a.  Grid  generation  and  WQM  and  HM  linkage. 

b.  WQM  and  HM  transport  comparisons  to  test  linkage. 

c.  Intratidal  and  intertidal  transport  tests.  ^ 

d.  WQM  calibration  for  DO  and  nutrient  simulation. 

e.  WQM  sensitivity  tests. 

/.  WQM  demonstration  of  use  for  evaluating  nutrient  load  reductions.  • 

The  DO  hypoxia  event  of  the  summer  of  1976  was  selected  for  model  applica¬ 
tion  since  this  period  had  relatively  abundant  data  and  was  a  period  of  high 
environmental  stress  (i.e.,  low  DO).  The  period  April  through  September  1976 
was  used  for  steps  a  through /above. 


Scope 


The  scope  of  die  study  reported  herein  was  restricted  to  studying  the  feasi- 
bility  of  modeling  this  complex  system.  Although  die  models  were  calibrated 
and  applied  for  various  sensitivity  conditions  and  nutrient  loading  scenarios, 
the  results  should  not  be  used  to  judge  die  effects  of  altered  nutrient  loadings 
on  eutrophication  and  water  quality.  This  study  did  not  include  sufficient 
detail,  accuracy,  and  degree  of  model  calibration/verification  required  for  a 

complete  nutrient  and  eutrophication  analysis.  • 

This  report  documents  die  application  steps  a  through  /  above.  The  steps 
were  conducted  concurrently  with  HM  and  WQM  development  Steps  a 
through  c  were  conducted  using  a  HM  grid  of  5  layers.  Step6  d  through  / 

were  conducted  using  an  enhanced  HM  grid  of  10  layers  with  minor  topogra-  q 

phic  changes. 

Chapter  2  gives  the  model  description,  and  Chapter  3  discusses  WQM  and 
HM  linkage  and  transport  comparisons.  Chapter  4  summarizes  input  data  and 
the  use  of  water  quality  observations  for  model  comparison.  Chapter  5  pre¬ 
sents  die  results  of  the  model-prototype  calibration,  sensitivity  tests,  and  • 

demonstration  scenarios.  Conclusions  and  recommendations  are  discussed  in 
Chapter  6. 
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2  Model  Description 


General 

The  NYB  HM  and  WQM  are  based  on  die  Chesapeake  Bay  modeling 
system  which  consists  of  three  interacting  models:  (1)  die  HM,  (2)  the  WQM, 
aid  (3)  a  bottom  sediment  diagenesis  model  (SM).  The  HM  is  an  improved 
version  of  CH3D  (Curvilinear  Hydrodynamics  in  Three  Dimensions)  developed 
by  Sheng  (1986)  for  the  U.S.  Army  Waterways  Experiment  Station  (WES). 

The  HM  was  extensively  modified  in  its  application  to  Chesapeake  Bay  (John¬ 
son  et  al.  1991a,b).  The  HM  operates  on  mi  intratidal  time  scale  (i.e.,  includes 
tidal  variations)  and  employs  a  curvilinear  or  boundary-fitted  planfonm  grid. 
Two  versions  of  CH3D  are  employed  by  WES:  one  version  has  stretched 
(sigma)  coordinates  for  die  vertical  dimension  while  die  other  has  fixed  verti¬ 
cal,  Cartesian  coordinates,  a  modification  for  the  Chesapeake  Bay  Study.  The 
sigma-stretched  version  was  used  for  the  NYB  study.  The  HM  was  exten¬ 
sively  tested  and  verified  during  the  Chesapeake  Bay  study  (Johnson  et  al. 
1991a). 

The  HM  and  WQM  are  operated  as  separate  modules.  Output  from  the 
HM  is  written  to  a  file  that  is  subsequently  used  as  input  by  the  WQM.  The 
WQM  can  accept  either  intratidal  or  intertidal  updates  for  hydrodynamics, 
which  are  converted  in  a  processor  built  into  the  HM  (Dortch  1990).  The 
procedure  is  computationally  efficient;  numerous  WQM  runs  can  be  executed 
without  recomputing  die  hydrodynamics. 

The  framework  of  the  WQM  (i.e.,  the  Chesapeake  Bay  WQM  and  the 
model  used  here)  is  based  on  the  finite  volume  approach,  similar  to  the 
U.S.  Environmental  Protection  Agency’s  (EPA)  WASP  (Water  Quality  Analy¬ 
sis  Simulation  Program)  model  (Ambrose,  Vandergrift,  and  Wool  1986). 

Finite  volume  transport  models,  also  referred  to  as  integrated  compartment  and 
multiple  box  models,  have  a  desirable  feature  in  that  they  can  be  linked  to 
finite  difference  and  finite  element  hydrodynamic  models. 

Significant  improvements  were  made  to  the  solution  schemes  of  the  Chesa¬ 
peake  Bay  model  compared  with  WASP  (Hall  and  Chapman  1982;  Ray  Chap¬ 
man  and  Associates  1988;  Hall  1990;  Cerco  and  Cole,  in  preparation).  The 
solution  recognizes  terms  in  the  horizontal  plane  and  vertical  direction.  Thus, 
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explicit  and  implicit  solution  schemes  are  used  in  the  horizontal  and  vertical 
dimensions,  respectively.  The  implicit  vertical  solution  can  greatly  reduce 
stability  constraints  on  time-steps.  A  third-order  accurate,  upstream-weighted 
differencing  scheme  (QUICKEST,  i.e.,  Quadradic  Upstream  Interpolation  for 
Convective  Kinematics  with  Estimated  Streaming  Terms,  Leonard  1979)  is 
used  for  the  horizontal  advective  terms  to  greatly  reduce  unwanted  numerical 
dampening.  Additionally,  the  WQM  time-step  is  variable  and  computed  within 
the  model  based  on  stability  requirements.  The  kinetic  routines,  which  include 
the  state  variables  and  their  interactions,  were  developed  during  the  Chesa¬ 
peake  Bay  Study  (Cerco  and  Cole,  in  preparation). 

The  WQM  and  SM  are  run  interactively  rather  than  coupled  indirectly  as 
the  HM  and  WQM.  The  SM,  described  by  DiToro  and  Fitzpatrick  (1992), 
simulates  decay  and  mineralization  (i.e.,  diagenesis)  of  organic  matter  depos¬ 
ited  in  the  sediments  and  the  exchange  of  nutrients,  DO,  and  other  substances 
between  the  sediments  and  overlying  water  column. 

This  study  included  the  application  of  the  HM  and  WQM,  but  the  SM 
module  was  not  used  because  little  or  no  data  required  by  the  SM  were  avail¬ 
able  for  this  region  when  the  study  was  initiated.  Adaptation  of  the  Chesa¬ 
peake  Bay  WQM  module  essentially  consisted  of  simulation  of  a  subset  of 
available  state  variables  with  minor  coding  changes  to  accommodate  exclusion 
of  the  SM  module  and  elimination  of  some  state  variables. 


Model  Adaptations 

The  Chesapeake  Bay  WQM  simulated  22  state  variables.  A  subset  of 
1 1  state  variables  were  selected  for  the  NYB  application.  Table  1  lists  the 
state  variables  simulated  in  the  NYB  study  and  the  naming  conventions  used  in 
the  NYB  and  Chesapeake  Bay  studies. 

Although  the  SM  module  was  not  used,  the  WQM  does  have  the  option  of 
specifying  benthic  fluxes  as  input  data.  Thus,  benthic  fluxes  for  sediment 
oxygen  demand  (SOD),  NH»-N,  and  N03-N  were  specified  in  this  study. 

Although  most  of  the  kinetic  rate  coefficients  were  specified  in  input  files, 
some  were  hardwired  in  die  computer  code.  Kinetic  coefficients  used  in  the 
NYB  study  are  tabulated  in  Appendix  A.  The  coefficients  labeled  “Hard¬ 
wired”  were  coded,  and  the  coefficients  labeled  “Input”  were  specified  in  input 
files. 
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3  Model  Linkage  and 
Transport  Comparisons 


The  linkage  and  testing  of  the  WQM  and  HM  modules  proceeded  through 
three  steps:  (1)  WQM  grid  generation  and  WQM  and  HM  grid  linkage, 

(2)  WQM  and  HM  transport  comparison,  and  (3)  intratidal  and  intertidal 
WQM  transport  tests.  The  linkage  and  testing  were  conducted  using  HM 
simulation  conditions  of  tidal  boundaries,  wind,  and  freshwater  flow. 


WQM  Grid  Generation  and  Linkage 

WQM  grid  generation  is  based  on  a  HM  grid  file  that  indudes  node  coor¬ 
dinates  and  cell  depths.  The  HM  grid  file  was  used  to  construct  a  map  file 
relating  die  HM  and  WQM  grid  cells.  The  WQM  grid  generation,  WQM  and 
HM  grid  linkage,  and  transport  tests  were  conducted  using  a  76  x  44  HM  grid 
consisting  of  3  layers.  Subsequent  modifications  of  the  HM  grid  consisted  of 
(1)  minor  changes  to  the  Connecticut  shoreline  and  Upper  Bay  and  (2)  the  use 
of  10  layers.  Plate  1  displays  the  initial  HM  grid  used  for  the  transport  tests 
while  Plate  2  displays  die  final  HM  grid  with  changes  to  die  Connecticut 
shoreline  and  Upper  Bay.  The  HM  grid  depicted  in  Plate  2  was  used  for 
WQM  calibration,  sensitivity  tests,  and  demonstration  scenarios. 

The  WQM  grid  was  a  direct  overlay  (i.e.,  one-to-one  cell  correspondence) 
of  the  76  x  44  cell  HM  grid.  However,  the  row  of  HM  cells  along  the  ocean 
boundary  were  not  included  in  the  WQM  grid  to  accomplish  proper  interfacing 
of  the  two  models.  Coding  appended  into  the  HM  used  die  mapping  informa¬ 
tion  to  initially  write  to  an  output  file  the  initial  cell  volumes  and  time- 
invariant  data  that  included  cell  surface  areas  and  horizontal  cell  lengths. 
During  the  HM  simulations,  temporally  averaged  flows  and  vertical  diffusiv- 
ities  were  computed  and  written,  along  with  cell  volumes  for  mass  continuity 
checks,  to  an  output  file.  The  details  of  these  procedures  are  explained  by 
Dortch  (1990). 

Proper  linkage  of  the  HM  and  WQM  was  ensured  through  mass  conserva¬ 
tion  tests.  Errors  in  the  linkage  are  easily  identified  by  observing  cell  mass 
deviations  that  cannot  be  accounted  for  through  elementary  mass  balance 


ChapWr  3  Modsl  Unfcaga  and  Transport  Comparisons 


calculations.  To  conduct  mass  conservation  tests,  output  data  from  the  HM 
were  used  to  simulate  three  conservative  tracer  scenarios:  (1)  initial  conditions 
(IQ  and  boundary  conditions  (BQ  of  tracer  concentration  set  to  zero;  (2)  IC 
and  BC  of  tracer  concentration  set  to  a  constant;  and  (3)  an  instantaneous  spot 
dump  of  tracer  mass  in  the  Apex  with  zero  IC  and  BC  for  tracer  concentration. 
In  the  first  test,  zero  concentrations  should  be  maintained  throughout  the  grid 
if  mass  is  conserved.  In  die  second  test,  a  constant  concentration  equal  to  the 
IC  and  BC  should  be  maintained  throughout  the  grid  if  mass  is  conserved.  In 
the  third  test,  no  mass  enters  or  leaves  the  grid,  so  die  sum  of  die  mass  within 
die  grid  should  equal  the  amount  of  mass  dumped.  Mass  conservation  was 
maintained  (within  the  accuracy  of  the  computer)  in  the  three  test  scenarios 
indicating  correct  linkage. 


HM  Versus  WQM  Transport  Comparison 

Transport  comparisons  between  the  HM  and  WQM  were  required  to  verify 
that  the  WQM  properly  represented  transport  provided  by  the  HM.  It  was 
assumed  that  die  HM  properly  simulated  transport  in  the  prototype.  The  salin¬ 
ity  state  variable  in  the  HM  and  WQM  was  used  to  simulate  a  conservative 
tracer.  Two  transport  comparisons  were  conducted:  (1)  continuous  tracer 
release  in  the  Hudson  River  and  (2)  continuous  tracer  release  in  die  Bight 


Hudson  River  continuous  release 

The  Hudson  River  experiment  consisted  of  maintaining  a  tracer  concentra¬ 
tion  of  10  units  in  the  column  of  cells  corresponding  to  the  junction  of  the 
Hudson  River  and  Upper  Bay.  Plate  3  displays  the  location  of  the  continuous 
release,  cells  sampled  for  time  history  plots,  and  cells  sampled  for  transect 
plots.  The  grid  displayed  in  Plate  3  is  a  window  of  the  HM  grid  shown  in 
Plate  1.  Only  a  window  of  the  grid  is  displayed  in  order  to  show  detail.  The 
solid  shaded  cell  represents  the  location  of  the  continuous  release,  the  open 
circles  represent  the  cells  sampled  for  time  series  [dots,  and  the  solid  lines 
represent  the  transects.  Plate  4  is  the  time  history  concentration  of  corre¬ 
sponding  HM  and  WQM  cells.  The  cells  displayed  are  surface  cells  located 
near  the  Brooklyn  Bridge  and  near  the  Verrazano  Narrows  Bridge  (WQM 
Cells  2422  and  2394,  respectively).  Hate  3  displays  the  two  transects  after 
S  days.  Transect  1  represents  the  surface  cells  extending  north  along  the  Hud¬ 
son  Canyon  into  Upper  Bay.  Transect  2  represents  a  transect  extending  from 
foe  Upper  Bay,  through  East  River,  and  into  Long  Island  Sound.  There  was 
"exact"  correspondence  between  HM  and  WQM  transport  using  equivalent 
time-steps,  flow  updates  at  time-step  intervals,  UPWIND  horizontal  advection, 
and  zero  vertical  advection.  The  HM  and  WQM  were  compared  using 
UPWIND  horizontal  advection  because  foe  HM  QUICKEST  code  required 
modification  in  order  to  handle  the  complex  boundaries  encountered  in  the 
Upper  Bay.  The  HM  QUICKEST  code  was  being  modified  during  the  trans¬ 
port  comparison. 
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Bight  continuous  ratoM* 

The  Bight  continuous  release  experiment  also  demonstrated  "exact"  corre¬ 
spondence  between  HM  and  WQM  transport  using  both  the  UPWIND  and 
QUICKEST  advecdon  schemes.  The  HM  and  WQM  were  compared  using  the 
QUICKEST  advectkm  scheme  because  HM  QUICKEST  boundary  formula¬ 
tions  did  not  affect  the  simulations.  The  continuous  release  occurred  in  the 
Bight  just  seaward  of  the  Rockaway  Point-Sandy  Hook  Transect  Plate  6 
shows  the  location  of  the  release  through  solid  shading  and  cells  sampled  for 
time  series  through  open  circles.  Plates  7  and  8  are  time  series  plots  at  corre¬ 
sponding  WQM  and  HM  cells  using  UPWIND  and  QUICKEST,  respectively. 
The  HM  and  WQM  transport  comparisons  indicate  that  the  WQM  duplicates 
HM  transport 


Intratidal  and  Intertidal  Transport  Comparison 

The  HM  flows  and  vertical  diffusivities  can  be  averaged  over  intratidal 
periods  (e.g.,  1  to  2  hr)  or  intertidal  periods  (i.e.  one  or  more  tidal  periods). 
The  intertidal  flows  are  averaged  in  a  manner  that  preserves  the  residual  cur¬ 
rents  associated  with  tidal  forcing.  The  procedure  involves  approximating  the 
Lagrangian  residual  velocities  as  the  sum  of  the  Eulerian  and  Stokes  velocities 
(Dortch  1990).  The  use  of  intertidal  hydrodynamics  reduces  time-step  restric¬ 
tions  in  the  WQM  and  greatly  reduces  the  amount  of  hydrodynamic  informa¬ 
tion  that  must  be  stored  and  read  in  during  WQM  execution,  thus  reducing 
WQM  execution  time  and  disk  storage  needs. 

During  the  sequence  of  mass  conservation  tests,  stable  time-step  sizes  for 
both  intratidal  and  intertidal  averaging  were  calculated.  On?  hour  intratidal 
averaging  permitted  a  WQM  time-step  of  17  min  and  a  Yl-'r  intertidal  aver¬ 
aging  permitted  a  52-min  time-step.  The  52-tnin  time-step  represents  a  lower 
bound  of  the  possible  time-step  because  the  calculated  time-step  in  the  WQM 
was  limited  to  a  maximum  of  60  min.  Intertidal  averaging  decreases  requited 
computer  execution  time.  However,  prior  to  the  use  of  intertidal  averaging,  it 
is  important  to  ensure  that  intertidal  averaging  duplicates  intratidal  transport. 

A  second  HM  and  WQM  Bight  tracer  comparison  was  conducted  using 
intratidal  (1-hr)  and  intertidal  (12.5-hr)  averaging.  Instead  of  a  continuous 
release,  a  spot  dump  was  simulated.  The  tracer  was  uniformly  injected  over 
the  depth.  Plate  9  displays  the  location  of  the  release  and  two  perpendicular 
transects  through  the  spot  dump.  Plate  10  represents  the  perpendicular  tran¬ 
sects  through  the  spot  dump  area  after  20  days  using  1-hr  intratidally  averaged 
flows.  Transect  1  represents  the  surface  cells  extending  north  along  foe  Hud¬ 
son  Canyon  info  Upper  Bay.  Transect  2  represents  a  transect  extending  from 
foe  New  Jersey  Coastline  to  Long  Island.  There  was  good  agreement  between 
HM  and  WQM  transport  using  1-hr  intratidal  averaging. 

Plate  1 1  demonstrates  the  same  perpendicular  transects  using  12.5-hr  inter¬ 
tidal  averaging.  Examination  of  Plate  11  reveals  that  a  phase  error  is 
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detectable  and  a  specification  of  S  m2  sec'1  horizontal  diffusion  in  the  WQM 
was  necessary  to  match  the  tracer  magnitudes.  The  tracer  tests  used  are 
extremely  dynamic;  the  tracer  front  is  being  advected  through  a  cell  in  a  few 
hours.  Any  time  averaging  awn  promises  the  transient  transport  information. 
Information  lost  resulted  in  the  appearance  of  phase  error  and  required  the 
addition  of  diffusion.  However,  even  considering  die  difficulty  of  simulating  a 
spot  tracer  dump  in  the  Bight,  the  tracer  experiment  confirmed  that  with  inter¬ 
tidal  averaging  the  WQM  could  duplicate  HM  transport  in  the  Bight  Except 
for  accidental  spills,  water  quality  constituents  in  the  Bight  do  not  exhibit  nidi 
strong  gradients  as  were  simulated  in  the  spot  dump.  Therefore,  intertidal 
averaging  is  appropriate  for  water  quality  simulation  in  the  Bight  For  simu¬ 
lating  spills  or  other  localized  conditions  with  extremely  sharp  gradients,  1-hr 
hydrodynamic  updates  should  be  used  in  the  WQM,  or  the  particle  tracking 
modd  that  was  developed  as  a  part  of  the  NYB  study  should  be  used. 

In  contrast  with  the  Bight  foe  current  dynamics  of  the  Upper  Bay  and  East 
River  display  greater  magnitude  and  nonlinearity.  A  series  of  experiments 
were  conducted  in  order  to  determine  if  intertidal  averaging  could  be  used  for 
the  Upper  Bay/East  River  and  select  foe  optimum  time-averaging  interval.  The 
experiments  consisted  of  using  1-hr  averaging  as  the  base  and  comping  dif¬ 
ferent  averaging  intervals.  The  comparisons  were  between  WQM  nd  did 
not  include  HM  tracer  comparisons. 

The  experiments  consisted  of  continuous  release  distributed  uniformly  from 
top  to  bottom  in  the  water  column  at  the  junction  of  foe  Hudson  River  and 
Upper  Bay.  Plate  12  displays  foe  location  of  foe  release  through  solid  shading 
and  cells  sampled  for  time  series  through  open  circles.  The  sampled  cells 
displayed  are  surface  cells  in  the  Upper  Bay  and  near  the  Verrazano  Narrow 
Bridge.  The  release  cell  is  at  foe  junction  of  foe  Hudson  River  and  Upper 
Bay.  Plates  13-17  represent  the  results  of  1-hr  intratidal,  4-hr  intratidal,  8-hr 
intratidal,  12-hr  intratidal,  and  12-hr  intertidal  averaging  for  foe  continuous 
release  experiment,  respectively.  Intratidal  is  used  here  to  refer  to  Eulerian 
averages,  and  intertidal  refers  to  Lagrangian  averages.  The  HM  was  executed 
for  90  days  with  flow  averaging  recorded  for  a  30-day  period  extending  from 
the  61st  through  90th  day  of  HM  execution.  The  WQM  simulated  a  300-day 
period  by  repeating  foe  30-day  HM  record.  Tracer  concentrations  were 
recorded  at  1-hr  intervals  for  plotting.  To  compare  the  different  flow¬ 
averaging  intervals,  the  plotting  data  were  Altered  with  a  running  12-hr 
average. 

Examination  of  die  plates  reveals  that  1-  and  4-hr  intratidal  averaging 
results  in  nearly  equivalent  asymptotic  results  while  8-  and  12-hr  intratidal 
averaging  results  in  a  decrease  in  magnitude.  The  12-hr  intertidal  averaging 
resulted  in  simulated  concentrations  that  asymptotically  approached  an  equilib¬ 
rium  value  slightly  exceeding  foe  1-  and  4-hr  values. 

Intratidal  transport  with  averaging  intervals  greater  than  4  hr  should  not  be 
used,  especially  in  the  bays  and  harbors  of  this  system.  Lagrangian  residual 
processing  should  be  used  for  intertidal  transport.  This  procedure  works  well 
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for  the  Bight  region  and  fairly  well  for  foe  bays  and  harbors  since  foe  intertidal 
transport  was  roughly  equivalent  to  foe  HM  and  foe  1-hr  transport  Although 
foe  intertidal  transport  through  the  transect  is  believed  to  be  properly  pre¬ 
served.  it  was  not  actually  tested.  In  retrospect,  it  should  be  verified  that  foe 
first-order,  Lagrangian  residual,  intertidal  transport  procedure  does  yield  the 
proper  net  flux  of  material  through  foe  transect  Such  tests  should  be  con¬ 
ducted  during  the  initial  phases  of  any  future  Bight  modeling. 
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4  Input  Data 


General 

The  eutrophication  WQM  required  six  types  of  input  data: 

a.  Hydrodynamic  transport  information. 

b.  Boundary  conditions  and  external  loadings  for  modeled  constituents, 

c.  Meteorological  data. 

d.  Initial  constituent  concentrations. 

e.  Kinetic  coefficients  for  transfers  and  transformations. 

/.  Observed  water  quality  data  for  calibration. 

Hydrodynamic  transport  information  was  stored  during  HM  simulation  as 
discussed  in  Chapter  3.  Boundary  conditions,  external  loadings,  initial  condi¬ 
tions,  and  water  quality  data  for  calibration  were  extracted  from  existing  water 
quality  observation  data  sets.  Meteorological  data  were  obtained  from  weather 
station  observations.  Kinetic  coefficients  were  adjusted  during  model  calibra¬ 
tion,  but  were  kept  within  expected  bounds  based  upon  previous  modeling 
experience  and  values  reported  in  the  literature. 

WQM  calibration,  sensitivity  tests,  and  demonstration  scenarios  were  con¬ 
ducted  with  the  grid  depicted  in  Plate  18,  which  included  10  vertical  layers. 
Plate  19  is  die  transformed  grid  corresponding  to  the  physical  grid  shown  in 
Plate  2.  The  transformed  grid,  which  is  used  in  the  HM  solution,  is  included 
to  more  easily  identify  the  location  of  a  cell  when  given  tire  HM  grid  coordi¬ 
nates.  The  WQM  grid  (Plate  18)  was  a  direct  overlay  of  die  76  x  44  cell  HM 
grid.  The  row  of  HM  cells  along  the  ocean  shelf  boundary  was  not  included 
in  the  WQM  grid  to  accomplish  proper  interfacing  of  the  two  models.  The 
WQM  used  254)10  computational  cells. 

Hydrodynamic  output  required  by  the  WQM  includes  initial  cell  volumes, 
surface  areas,  and  horizontal  cell  dimensions;  three-dimensional  flows  for  all 


12 


ChapfcrS  Input 


cell  frees;  and  vertical  eddy  diffusivities.  Time  of  each  hydrodynamic  update 
and  cell  volumes  for  each  update  are  also  output  for  continuity  checks. 

Water  quality  observations  were  used  by  the  WQM  to  specify  initial  and 
boundary  conditions  and  external  loadings,  and  were  used  for  calibration  com¬ 
parisons.  Readily  available  water  quality  data  sets  collected  during  1975-1976 
within  die  NYB  were  compiled,  documented,  and  furnished  in  digital  form  by 
Creative  Enterprises,  Inc.,  and  Han  A  Associates,  Inc.  (1990).  The  use  of 
water  quality  data  is  described  in  more  detail  in  the  next  section.  Kinetic 
coefficients  are  discussed  in  Chapter  5. 

Daily  averaged  solar  radiation,  water  equilibrium  temperatures,  and  coeffi¬ 
cients  of  surface  beta  exchange  were  computed  from  meteorological  data. 

Solar  radiation  values  were  used  in  the  algal  growth  functions,  and  the  equi¬ 
librium  temperatures  and  die  surface  heat  exchange  coefficients  were  used  as 
input  for  temperature  simulation  (Edinger,  Brady,  and  Geyer  1974).  The  com¬ 
putations  of  these  three  variables  from  meteorological  data  were  based  on  a 
Heat  Exchange  Program  (Eiker  1977).  The  U.S.  Air  Force  Environmental 
Technical  Applications  Center  provided  meteorological  data  in  the  form  of 
Tape  Deck  1440  WBAN  Hourly  Surface  Observations  for  JFK  Airport.  Algal 
growth  functions  also  required  estimates  of  photoperiod.  Photoperiod,  the 
daylight  fraction  of  a  24-hr  day,  was  computed  using  a  formula  in  Stoddard 
(1983). 


Use  of  Water  Quality  Observations 

The  use  of  observed  water  quality  data  for  initial  conditions  and  boundary 
conditions  is  described  in  this  section.  The  special  treatment  of  observations 
for  WQM  calibration  comparisons  is  described  in  Chapter  5  (Application 
Results). 


Initial  conditions 

The  WQM  simulations  extended  from  15  Apr  through  30  Sep  1976.  Initial 
conditions  were  specified  for  each  of  the  WQM  state  variables  and  for  each 
cell  at  the  initiation  of  die  simulation.  The  water  quality  data  collected  during 
April  1976  were  used  to  specify  initial  conditions.  Because  of  the  sparsity  of 
data,  initial  conditions  had  to  be  specified  on  a  regional  basis.  This  was 
accomplished  with  an  overlay  grid. 

A  coarse  overlay  grid  was  superimposed  on  the  transformed  WQM  grid 
(Plate  20).  The  15  coarse  grid  cells  roughly  correspond  to  available  data  and 
the  study  area’s  bathymetry.  Table  2  lists  the  number  of  surface  cells  included 
in  the  coarser  grid  cells  and  the  average  depth  of  the  coarser  grid  cells. 
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34 

12 

3 

236 

22 

4 

140 

18 

S 

186 

24 

6 

277 

27 

7 
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43 
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130 

38 
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42 

10 

168 

53 

11 

264 

70 

12 

182 

132 

13 

120 

118 

14 

156 

124  1 

I  15 

156 

143  | 

The  coarse  overlay  grid  also  had  10  vertical  layers.  The  observed  data 
were  assigned  to  each  of  the  ISO  cells  (IS  surface  coarse  grid  cells  x  10  lay¬ 
ers)  based  on  its  horizontal  location  and  depth  of  collection.  The  observed 

data  for  each  coarse  grid  cell/layer  were  averaged  to  obtain  a  single  value  for  • 

that  coarse  cell/layer.  Measured  data  were  not  available  for  assignment  to  each 
of  the  ISO  cells;  however,  each  coarse  cell  required  an  assignment  for  each 
state  variable.  The  following  procedure  was  used  to  assign  missing  values. 

Examination  of  observed  data  revealed  two  vertical  patterns  of  distribution 

which  were  used  for  vertical  assignments.  One  pattern,  exemplified  by  all  the  q 

state  variables  except  plankton,  varied  with  depth  from  a  surface  value  to  a 

value  that  remained  constant  through  most  of  the  water  column.  The  depth  at 

which  values  approached  constancy  corresponded  to  the  depth  of  die  pycno- 

cline.  The  second  pattern  was  characteristic  of  the  phytoplankton.  Maximum 

values  occurred  near  pycnocline  depth  with  minimum  values  at  the  bottom. 

Based  on  these  observed  vertical  distributions,  missing  horizontal  cells  for  the  • 

surface,  middle,  and  bottom  layers  for  the  plankton  and  surface  and  bottom 
layers  for  the  other  state  variables  were  assigned  to  the  coarse  grid.  Assigned 
values  were  selected  to  maintain  existing  gradients.  Values  were  interpolated 
for  layers  that  did  not  have  data  as  explained  in  the  next  paragraph.  The  ini¬ 
tial  condition  values  used  for  the  ISO  coarse  overlay  cells  are  tabulated  in  ^ 

Appendix  B.  Observed  values  are  underlined.  The  same  overlay  values  were 
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then  used  for  die  initial  conditions  of  each  computational  cell  within  the  coarse 
overlay  cell  and  layer. 

The  initial  condition  file  listed  in  Appendix  B  specifies  three  values  for  net 
plankton  and  nanoplankton  corresponding  to  surface,  pycnodine,  and  bottom. 
All  the  other  constituents  specify  two  values,  representing  surface  and  bottom. 
For  constituents  with  two  values,  the  above  pycnodine  layers  were  assigned 
values  varying  linearly  between  the  surface  and  bottom  value.  The  below- 
pycnocline  layers  were  assigned  the  bottom  value.  For  plankton,  the  layer 
assignments  were  assigned  through  a  parabolic  function.  For  layers  above  the 
pycnodine,  the  pycnodine  value  was  reflected  through  the  surface;  for  layers 
below  the  pycnodine,  the  pycnodine  value  was  reflected  through  tire  bottom. 

Reflection  needs  some  explanation.  Three  noncoincident  points  uniquely 
specify  tire  coefficients  of  a  parabola.  The  reflection  of  the  pycnodine  value 
through  tire  bottom  implies  that  an  imaginary  point  of  value  equal  to  the  pyc¬ 
nodine  value  is  spedfied  below  the  bottom.  The  distance  below  the  bottom  of 
the  imaginary  point  is  equal  to  the  distance  between  the  pycnodine  and  bot¬ 
tom.  A  parabola  is  then  fitted  using  the  pycnodine  value,  bottom  value,  and 
the  imaginary  value  located  bdow  the  bottom.  Reflection  through  the  surface 
is  similar,  the  three  points  represent  the  pycnodine  value,  surface  value,  and 
an  imaginary  value  located  above  the  surface.  The  above-surface,  imaginary 
point  has  a  value  equal  to  tire  pycnodine  value  and  is  located  above  the  sur¬ 
face  a  distance  equal  to  the  distance  between  the  pycnodine  and  surface.  The 
reflection  procedure  was  used  to  ensure  maximum  phytoplankton  at  pycnodine 
depth  with  rapid  attenuation  to  surface  and  bottom  values.  The  pycnodine 
depth  was  uniformly  specified  as  20  m  for  cell  column. 


Boundary  conditions 

Boundary  condition  files  include  the  temporal  and  spatial  specification  of 
constituent  concentrations  for  the  seaward  boundary  and  Hudson  River,  and 
specification  of  external  mass  loadings  for  all  other  sources  of  nutrients/water 
quality.  Boundary  condition  files  are  also  used  to  specify  meteorological  data 
and  sediment-water  fluxes.  Monthly  updates  were  used  for  tire  seaward  boun¬ 
dary  and  the  Hudson  River,  and  the  external  mass  loadings  were  held  constant 
These  update  frequencies  were  considered  adequate  for  the  purpose  of  this 
study,  that  is,  a  feasibility  study.  For  a  comprehensive  eutrophication/nutrient 
analysis,  more  frequent  updates  should  be  used,  such  as  daily  or  weekly 
updates  for  the  Hudson  River  and  monthly  updates  for  external  loadings. 

Seaward  boundary.  Measured  data  collected  seaward  of  the  shelf  break 
boundary,  seaward  of  Atlantic  Gty,  and  seaward  of  Montauk  Point  were  used 
to  prepare  seaward  boundary  conditions.  The  seaward  boundary  was  divided 
into  seven  segments,  shown  in  Plate  21.  Segment  8  refers  to  the  Hudson  River 
boundary  segment  which  is  described  below.  The  boundary  condition  values 
are  tabulated  in  Appendix  C.  Observed  values  are  underlined.  The  seaward 
boundary  condition  values  were  specified  similarly  to  the  initial  condition  files; 
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(1)  assign  observed  data  to  each  of  (he  840  space-time  slots  (7  surface  seg¬ 
ments  x  10  layers  x  12  months)  and  (2)  assign  layer  values  based  on  surface, 
pycnocline,  and  bottom  values  as  described  above.  Then  each  computational 
cell  is  assigned  die  space-time  slot  value  based  on  the  simulation  clock  time 
and  the  cell  location  with  respect  to  die  segments  Kid  layers.  Each  entry  in 
Appendix  C  corresponds  to  a  surface,  pycnocline,  and  bottom  space -time  slot 
Examination  of  Appendix  C  reveals  the  sparsity  of  data  available  to  fill  the 
space-time  slots.  Empty  space-time  slots  were  filled.  The  criterion  used  for 
filling  empty  space-time  slots  was  maintenance  of  spatial  and  temporal 
continuity. 

Hudson  River  upstream  concentrations  and  loadings.  Water  quality 
constituents  can  enter  the  Hudson  River  as  "upstream"  inputs  (i.e.,  upper 
watershed  runoff)  and  local  point  and  nonpoint  source  loads  (i.e.,  municipal 
and  industrial  discharges  and  tributaries  that  enter  the  Hudson  downstream  of 
the  upstream  model  limits  for  the  river).  The  Hudson  River  upstream  inputs 
were  specified  in  two  ways,  as  concentrations  and  mass  loadings.  Tempera¬ 
ture,  salinity,  and  DO  were  specified  as  monthly  varying  boundary  concentra¬ 
tions.  Temperature  values  were  calculated  monthly  averaged  equilibrium 
temperatures.  Dissolved  oxygen  values  represented  80  percent  saturation  at  the 
specified  temperature  and  salinity.  Salinity  was  specified  constant  at  20  parts 
per  thousand  (ppt).  The  Hudson  River  boundary  segment  is  noted  as  Seg¬ 
ment  8  on  Plate  21,  and  the  monthly  temperature,  salinity,  and  DO  concentra¬ 
tions  are  tabulated  in  Appendix  C. 

Due  to  data  limitations,  it  was  necessary  to  use  mass  loadings  for  the  other 
water  quality  constituents  entering  through  the  upstream  model  limits  for  the 
Hudson  River.  The  USGS  published  daily  measured  flows  at  Green  Island, 
downstream  from  the  Troy  and  Mohawk  River  locks  and  dams.  Mueller  et  al. 
(1976)  published  water  quality  data  collected  near  Poughkeepsie,  which  is 
downstream  from  Green  Island.  Mueller  used  drainage  basin  areas  and  their 
ratio  to  compute  a  factor  of  1.45  for  scaling  measured  flows  from  Green  Island 
to  estimated  flows  at  Poughkeepsie.  Measured  USGS  flows  at  Green  Island, 
Meuler’s  measured  water  quality  values  (Table  3),  and  the  scale  factor  of  1.45 
were  used  to  compute  daily  Hudson  River  tributary  mass  loads. 

Simulation  period  and  annual  averages  of  the  computed  daily  loads  for 
1976  are  compared  in  Table  4  to  annual  average  daily  loads  estimated  by 
Mueller  et  al.  (1976).  Table  4  reveals  that  the  computed  loads  used  in  this 
study  are  nearly  equivalent  to  those  reported  by  Mueller  et  al.  (1976).  More 
accurate  procedures  for  estimating  Hudson  River  loads  should  be  explored  for 
a  comprehensive  eutrophication/nutrient  analysis. 

Constant  external  constituent  loads.  Constant  external  constituent  loads 
(kilograms  per  day),  arising  from  local  runoff,  point  source  discharges,  and 
dump  sites,  were  defined  as  follows: 
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New  Yoik-New  Jersey  Harbor  and  estuaries 
(referred  to  as  Transect  region) 

New  Jersey  (NJ)  Coast 
Long  Island  (LI)  Coast 

Sewage  sludge  dump  site  • 

Dredged  material  dump  site 

The  constant  external  constituent  loads  used  in  this  study  are  listed  in  Table  5 
and  were  extracted  from  Mueller  et  al.  (1976).  The  Transect  external  loads 
included  wastewater  discharges  and  Mueller's  ungaged  runoff  loads  (Hudson 

River  drainage  below  Poughkeepsie,  New  York  City  area,  and  New  Jersey  ® 

ungaged  area  extending  north  from  Sandy  Hook  along  Raritan  Bay  to  the 

Hackensack  River).  The  New  Jersey  Coast  external  loads  included  wastewater 

discharges  and  surface  runoff.  The  Long  Island  Coast  external  loads  included 

wastewater  discharges,  surface  runoff,  and  groundwater  mass  loads. 

The  constant  external  loads  were  inserted  uniformly  throughout  foe  * 

10  layers.  The  Transect,  NJ  Coast,  and  LI  Coast  loads  were  distributed  uni¬ 
formly  in  cells  adjacent  to  the  shoreline.  Sewage  sludge  dump  site  loads  were 
inserted  in  foe  cells  corresponding  to  44°22'30‘-40°2500'  N  Latitude, 

73o4T30"-73°45’00‘  W  Longitude.  The  dredged  material  dump  site  loads  were 

inserted  in  foe  cells  corresponding  to  40°23'38"-40°2T28"  N  Latitude,  • 

73°5 1  ’28"-73°50’00"  W  Longitude. 

Sediment/water  column  boundary.  The  sediment/water  column  interface 
boundary  condition  was  specified  as  temporally  and  spatially  constant  for  SOD 

and  ammonium  (NH«-N)  release.  Both  SOD  and  NH«-N  releases  were  modu-  ^ 

lated  through  temperature  and  DO  concentration  of  overlying  water.  The 
sediment  ammonium  retease  at  reference  temperature  and  with  overlying  oxy¬ 
genated  water  was  0.043  g  N  m'2  day1.  Ammonium  release  increased  with  a 
decrease  in  DO.  The  SOD  and  ammonium  release  functions  are  described  in 
Appendix  A. 
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An  examination  of  available  SOD  measurements  for  both  1974  and  1975 
revealed  that  data  were  spatially  limited  to  collection  in  the  Apex  only 
(Plate  22).  The  reported  SOD  data  had  been  normalized  to  20  °C.  The  SOD 
data  woe  further  normalized  to  nonoxygen  limited  values  by  assuming  a  DO 
half-saturation  constant  of  2.0  g  O,  m  ’  and  using  recorded  overlying  water  DO 
concentration  values.  Table  6  displays  the  measured  SOD  values  and  the 
number  of  observations.  The  rows  and  columns  represent  the  corresponding 
HM  (i  j)  cells.  The  arithmetic  average  of  all  normalized  SOD  measurements 
was  0.88  g  Qj  m4  day'1.  Because  of  the  limited  spatial  extent  of  SOD  mea¬ 
surements  and  the  feasibility  aspect  of  this  study,  a  spatially  and  temporally 
constant  SOD  value  of  1.0  was  used.  A  sensitivity  test  of  SOD  is  described 
later  in  this  report 
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5  Application  Results 


The  summer  hypoxia  event  of  1976  was  selected  for  the  WQM  application. 
The  1976  hypoxia  event  was  selected  because  Dr.  Andrew  Stoddard  (Stoddard 
1983)  conducted  a  comprehensive  and  thorough  search  and  analysis  of  existing 
data  and  literature  concerning  die  eutrophication  and  oxygen  depletion  in  the 
New  York  Bight,  and  developed  and  applied  a  mathematical  model  of  oxygen 
depletion  of  the  1976  hypoxic  event  in  the  New  York  Bight-  Dr.  Stoddard’s 
effort  provided  a  basis  for  the  evaluation  of  die  WQM  application.  The  model 
was  first  calibrated  (i.e.,  parameters  were  adjusted)  against  observed  data  for 
dus  period.  Sensitivity  studies  were  then  conducted  with  the  model  to  investi¬ 
gate  the  importance  of  some  parameters  and  boundary  conditions.  The  WQM 
was  run  with  variations  in  nutrient  loadings  to  demonstrate  model  utility  and  to 
obtain  a  better  sense  of  the  effect  of  loads  on  die  system. 

The  WQM  was  also  used  to  investigate  die  causes  of  the  New  Jersey  near¬ 
coast  hypoxia,  die  event  that  occurred  during  die  1976  hypoxia.  Finally,  LIS 
external  loads  were  included  and  results  were  re-evaluated.  Water  quality  data 
for  LIS  external  constituent  loads  were  not  readily  available  at  die  initiation  of 
the  study,  and  observations  for  calibration  within  US  were  not  provided  in  die 
observation  data  sets.  Calibration,  sensitivity  tests,  and  load  reduction  scenar¬ 
ios  woe  conducted  without  LIS  external  loads.  An  investigation  was  ccm- 
ducted  to  examine  die  impact  of  LIS  external  loads  on  the  DO  resources  of  the 
Bight 

The  HM  simulation  extended  from  1  Apr  through  30  Sep  1976,  for 
183  days.  The  HM  output  was  averaged  over  12.5-hr  intervals  beginning 
13  Apr  and  written  to  a  file  for  subsequent  input  to  die  WQM.  The  183-day 
HM  simulation  required  31  hr  of  computation  time  on  a  CRAY  Y-MP.  The 
168-day  (15  Apr  through  30  Sep)  WQM  simulation  required  2  hr. 


Model-Prototype  Calibration 

The  temporal  and  spatial  sparsity  of  prototype  observations  for  water  qual¬ 
ity  made  model-prototype  comparisons  difficult  Monthly  temporal  averaging 
of  all  prototype  measurements  and  model  simulations  was  done  to  mitigate  the 


temporal  sparsity  of  prototype  observations.  Several  additional  techniques 
were  used  to  provide  meaningful  model-prototype  comparisons: 


II.  Scatter  plots 

2.  Regional  and  point  comparisons 
a.  Spatial  average 

!  (1)  Time  series  (dots 

|  (2)  Depth  profiles 

b.  Depth  profiles  at  specific  locations  (point  comparisons) 

*  3.  Transect  {dots 

4.  Statistics 

An  explanation  of  each  of  die  above  comparison  methods  is  discussed  in  its 
|  respective  section  along  with  die  results.  The  model-prototype  calibration 

(dots  are  displayed  in  Appendix  D. 
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Scatter  plots 

Measured  data  were  aggregated  by  region  (Apex,  NJ  Coastal  Zone,  etc.), 
layer,  and  month.  Each  measured  datum  was  associated  with  die  correspond¬ 
ing  aggregated  simulated  result  and  displayed  by  plotting.  The  aggregation 
procedure  is  discussed  more  completely  in  the  next  section.  The  scatter  (dots 
display  predicted  versus  observed  data  for  all  regions  and  were  effective  for 
gross  qualitative  calibration 

Plate  D1  displays  the  scatter  plots  for  the  calibration  run.  Examination  of 
Plate  D1  provides  a  quick  reference  as  to  ow  well  the  model  matches  obser¬ 
vations.  Plate  D1  indicates  that  simulated  *t  plankton,  DOC,  and  POC  were 
generally  underestimated.  However,  Appendix  B,  Initial  Conditions  Constitu¬ 
ent  Concentrations,  and  Appendix  C,  Boundary  Constituent  Concentrations, 
indicate  that  few  net  plankton,  DOC,  and  POC  measurements  were  available 
fra  initial  condition  and  boundary  specification.  Future  studies  should  attempt 
to  fill  this  data  gap. 


Regional  comparison 

There  were  few  stations  that  were  frequently  sampled  at  the  same  location. 
One  technique  used  to  provide  meaningful  model-prototype  comparisons  was 
spatial  aggregation  to  provide  larger  sample  sizes.  This  technique  consisted  of 
computing  monthly  averages  of  multiple  station  values  that  were  contained 
within  a  region.  The  NYB  was  subdivided  into  1 1  regions  as  shown  in 
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Plate  23.  The  segmentation  scheme  is  a  modification  but  reflects  the  segmen¬ 
tation  used  by  NOAA  (1979).  No  data  were  provided  (by  Creative  Enterprises 

and  Han  A  Associates  1990)  for  US,  rod  data  obtained  in  die  eastern  • 

(upcoast)  regions  were  used  for  upcoast  boundary  conditions.  Therefore,  die 

eastern  regions  were  not  included  in  comparisons  of  model  versus  observed 

data,  bid  US  results  (Region  11)  were  included  in  the  sensitivity  studies  to 

observe  the  test  effects  on  the  Sound.  Values  for  the  number  of  surface  cells 

and  average  depth  of  the  11  spatial  aggregates  are  given  in  Table  7. 


Table  7 

Regional  Spatial  Aggregatee 

nepoo 

No.  Surtaoa  Ca8a 

Avaraga  Dapth,  n 

1.  Apax 

184 

22. 

2.  NJ  nearshore  (<30  m) 

213 

21. 

3.  NJ  micfehatf  (30-60  m) 

262 

42. 

4.  NJ  shaN  <>60  m) 

123 

08. 

5.  HR  mkfchsH  (30-00  m) 

110 

48. 

6.  HR  shaft  (>60  m) 

00 

152. 

7.  LI  nearshore  (<30  m) 

80 

27. 

8.  LI  mktehaH  (3060  m) 

382 

46  | 

|  0.  Ushatf(>60m) 

313 

134.  | 

J  10.  Transact 

77 

6.  1 

11.  U  Sound 

170 

17.  j 

Note:  NJ  is  Now  Jersey.  HR  is  Hudson  River,  and  LI  is  Long  Wand.  | 

The  use  of  regional  aggregates  can  also  result  in  biases  that  confound 
model  calibration.  The  model  results  were  averaged  for  all  model  time-steps 
within  a  given  month  and  for  all  model  cells  within  the  region,  thus  producing 

a  true  monthly,  regional  average.  However,  the  prototype  data  are  far  from  q 

representing  true  monthly,  regional  averages  because  of  the  sparsity  of  data. 

Therefore,  specific  sampling  locations  were  selected  to  compare  measured  and 
simulated  data.  The  sample  locations  selected  for  point  comparisons  were 
based  on  the  availability  of  data  (e.g.,  temperature,  salinity,  and  dissolved 
oxygen).  Plate  23  displays  the  relationship  between  the  spatial  regions  and 

selected  point  comparisons.  The  regions  are  outlined  with  bold  lines  and  the  • 

selected  points  marked  with  solid  circles.  Table  8  lists  the  corresponding 
surface  HM  and  WQM  cell  and  depth  of  the  selected  points. 

The  ordered  (i  j)  following  the  regional  name  in  Table  8  represents  the 
corresponding  HM  surface  cell.  The  regional  spatial  comparisons  were  dis-  ^ 

played  through  depth  profiles  for  various  months  and  time  series  for  a  layer 
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Tibfo8 

8d9did  Point  Compartoono 

NH  Cal 

WON  Cal 

Dap*,  at 

1.  Apa* 

(31.32) 

2003 

21. 

2.  NJ  na— how  (<80  m) 

<  4.10) 

1201 

22. 

S.  NJ  ihWMmM  (SO-SO  m) 

(22,10) 

1067 

44. 

4.  NJ  thatf  (>S0  m) 

(17.  4) 

104 

104. 

5.  HR  mfcfcha*  (3040  m) 

(30.20) 

1361 

42. 

6.  HR  shaV  (>60  m) 

(28,11) 

80S 

as. 

7.  U  naarahora  (<30  m) 

(48,20) 

1746 

34. 

8.  U  midahatf  (30-60  m) 

(40.20) 

1377 

36. 

S.  U  ahai  (»80  m) 

(43,  0) 

436 

73. 

10.  Tranaact 

(31,30) 

2300 

0. 

(i.e.,  near  surface  or  near  bottom).  The  point  comparisons  were  displayed 
through  depth  profiles  at  points  corresponding  to  the  spatial  comparisons. 

The  depth  profiles  display  simulated  and  measured  constituent  values  as  a 
function  of  model  layer.  It  is  important  to  note  for  plot  interpretation  that  the 
Tables  7  and  8  list  the  depths  for  each  region  and  point  that  can  be  used  to 
compute  depths.  For  example,  the  N)  Nearshore  region  has  an  average  depth 
of  21  m  while  the  NJ  Shelf  region  has  an  average  depth  of  %  m  (Table  7). 
The  WQM  vertical  grid  has  10  equally  thick  layers;  therefore,  the  center  of  die 
bottom  layer  of  die  NJ  Nearshore  represents  a  depth  of  20  m  while  the  center 
of  the  bottom  layer  of  the  NJ  Shelf  represents  a  depth  of  91  m.  The  measured 
prototype  data  are  presented  at  the  appropriate  layer  depths.  The  results  of  the 
depth  profiles  follow  the  time  series  plots. 

lime  aeries.  The  time  series  plots  represent  the  monthly  averaged  surface 
and  bottom  layer  model  and  prototype  data  for  each  region  between  April  and 
September  1976.  Plates  D2  through  D9  display  regional  time  series  for  con¬ 
stituents  as  follows: 
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Plate  D2  indicates  that  simulated  temperatures  approximate  measured  val¬ 
ues.  Time  series  of  simulated  and  measured  DO  (Plate  D9)  also  indicate  rela¬ 
tively  close  agreement  However,  die  surface  and  bottom  layer  simulated  DO 
exceeded  measured  values  in  the  Transect  region.  Furthermore,  simulated  bot¬ 
tom  DO  in  the  Transect  region  frequently  exceeded  surface  values.  The  exces-  • 

rive  DO  computed  in  the  Transect  region  is  attributed  to  the  nanoplankton 
blooms  tint  occurred  in  the  modeled  Transect  region.  Nanoplankton  photosyn¬ 
thesis,  with  resulting  DO  production,  was  high  in  Raritan  Bay  and  occurred 
throughout  the  model  depth.  The  column  depths  of  the  cells  in  Raritan  Bay 
were  3  to  6  m.  It  is  not  known  how  excessive  the  model  algal  production  is 

since  there  are  no  nanoplankton  observations  in  die  Transect  region.  It  is  ® 

speculated  that  model  production  should  be  decreased  in  this  region  through 

light  attenuation  associated  with  suspended  solids.  It  is  feasible  for  surface 

DO  concentrations  to  be  lower  than  bottom  values  since  the  surface  layer 

would  rapidly  equilibrate  toward  saturation  through  release  to  the  atmosphere. 

A  reasonable  hypothesis  is  that  prototype  bottom  layer  DO  concentration  was  • 

less  than  simulated.  Suspended  solids  would  attenuate  light  transmission, 
resulting  in  little  photosynthetic  DO  production  in  the  lower  layers. 

Plates  D5  and  D6  (DOC  and  POC  time  series)  indicate  temporally  sparse 
measured  data.  The  large  average  values  of  POC  in  the  transect  are  due  to  die 

large  simulated  nanoplankton  concentrations  (Plate  D4).  Algal  carbon  is  added  * 

to  the  state  variable  POC  for  display  since  die  sum  of  the  state  variables  POC, 

net  plankton  cartoon,  and  nanoplankton  carbon  approximates  die  laboratory 

measurement  of  POC.  Net  plankton  concentrations  were  low  and  visually 

appear  diminutive  (Plate  D3)  because  of  the  selection  of  die  ordinate  values  for 

display;  die  maximum  value  of  20  mg  Chi  m  3  was  selected  because  of  die  one  • 

measurement  observed  in  May  in  the  bottom  layer  of  die  Hudson  River  Mid¬ 
shelf  region. 

The  WQM  underpredicted  DOC  in  the  Apex  yet  overpredicted  DOC  in  the 
Transect  The  excessive  DOC  in  the  Transect  is  the  result  of  hydrolysis  of  die 

large  values  of  POC  computed,  due  to  nanoplankton  computed.  However,  the  ® 

model  underpredicts  DOC  Bight-wide.  The  underprediction  is  reflected  in  the 


Chapter  S  Appfcaton  Fteauhi 


25 


DOC  scatter  plot  (Plate  Dl).  It  is  believed  Oat  the  ocean  contains  refractory 
DOC  that  degrades  very  slowly.  The  model  has  only  one  DOC  variable.  To 
more  accurately  simulate  DOC,  the  model  should  have  two  DOC  variables, 
labile  and  refractory.  Labile  DOC  degrades  much  faster  than  refractory  DOC. 
Presently,  the  model  DOC  degrades  at  rates  more  representative  of  labile  DOC. 

Examination  of  Plate  D7  indicates  a  spike  in  ammonia  during  May  in  the 
Transect.  There  were  no  measured  ammonia  values  in  the  Transect  during  the 
period,  but  examination  of  Plate  D8  reveals  a  spike  in  measured  nitrate  during 
May.  which  tends  some  credibility  to  an  ammonia  spike  since  ammonia  oxi¬ 
dizes  to  nitrate. 

Profiles.  Plates  DIO  and  Dl  1  display  regionally  averaged,  vertical  profiles 
for  temperature  and  DO,  respectively,  for  the  month  of  August  In  the  shallow 
waters,  such  as  the  Transect  region,  both  measured  and  observed  temperature 
and  DO  indicate  a  well-mixed  water  column  In  the  deeper  waters,  the  vertical 
profiles  of  temperature  indicate  similar  surface  and  bottom  temperatures  for 
both  measured  and  simulated  data,  but  the  depths  of  the  thennocline  are  gener¬ 
ally  deeper  for  the  measured  data  (e.g.,  HR  Shelf  and  NJ  Shelf).  It  is  hypothe¬ 
sized  that  these  disagreements  are  the  result  of  too  little  vertical  diffusion  in 
the  surface  layers  in  the  deeper  waters.  The  vertical  diffusion  used  in  the 
WQM  is  an  output  variable  from  the  HM.  Earlier  HM  output  had  produced 
too  much  vertical  mixing,  whereas  these  final  results  seem  to  produce  too  little 
mixing.  Clearly,  future  efforts  should  focus  on  more  accurate  representation  of 
vertical  diffusion  in  the  HM  through  one  of  foe  higher  order  turbulence  closure 
models  that  are  available. 

The  discrepancies  between  measured  and  simulated  vertical  profiles  of 
temperature  and  DO  were  similar  (Plates  DIO  and  Dl  1).  Notably,  foe  simu¬ 
lated  DO  was  about  a  factor  of  two  greater  than  that  measured  in  foe  Transect, 
due  to  excessive  algal  production  in  this  region.  The  reason  for  this  problem 
was  discussed  above  in  foe  explanation  of  foe  DO  time  series  plots  (Plate  D9). 
Also,  as  stated  earlier,  foe  sparsity  of  DO  observations  in  foe  Transect  region 
may  make  comparisons  of  regionally  averaged  DO  meaningless  for  this  region. 

Vertical  profiles  in  August  1976  were  compared  to  simulated  results  at 
specific  stations,  rather  than  using  regional  averages,  as  shown  in  Plates  D12 
and  D13  for  temperature  and  DO,  respectively.  Comparison  of  foe  computed 
area  and  point  profiles  indicates  similar  vertical  distributions.  However,  foe 
computed  thennocline  depth  was  deeper  at  foe  point  (28,11),  HR  Shelf,  than 
foe  average  over  foe  region,  thus  comparing  more  closely  with  foe  average 
measured  profile  for  this  region  The  simulated  bottom  DO  values  for  point 
comparisons  were  less  than  foe  area  averages. 


Transact! 

Transect  {Hots  are  another  procedure  for  comparing  model-prototype  results. 
Plate  24  shows  the  five  transects  that  were  used.  Their  location  and  width 


26 


Chapter  S  AppfcaSon  Oeauhi 


were  rejected  to  maximize  the  amount  of  available  prototype  data.  With  this 
technique,  model  results  and  observations  were  avenged  over  each  month  and 
across  the  transect  The  results  were  plotted  m  monthly  avenged  surface  and 
bottom  layer  concentration  versus  distance  from  the  seaward  boundary  along 
the  transect 

Plates  D14  and  D1S  display  surface  and  bottom  transect  {dots  for  August 
1976  for  temperature  and  DO,  respectively.  The  temperature  transects  indicate 
that  simulated  temperatures  match  well  the  measured  temperatures  for  August 

Examination  of  DO  transects  (Plate  D15)  reveals  that,  overall,  the  computed 
and  observed  DO  match  fairly  well  except  for  foe  bottom  layer,  nearshore  of 
Transect  1,  where  computed  values  exceeded  measured  values.  Several  things 
could  have  led  to  this  discrepancy,  such  as  monthly  averaging  of  model  results 
and  the  use  of  constant  shoreline  and  point  source  loads  and  monthly  updates 
for  boundary  loads.  Examination  of  foe  submonthly,  time-varying  simulated 
DO  (at  model  time-step  intervals)  for  foe  bottom  layer  of  Transect  1  indicated 
that  simulated  DO  was  less  than  1.0  g  m5  at  times.  However,  foe  model 
results  in  Plate  D1S  represent  monthly  averages.  Monthly  averaging  may 
provide  adequate  resolution  for  the  detection  of  seasonal  trends,  but  it  attenu¬ 
ates  extremes  that  occur  over  a  few  days  aid  may  foil  to  detect  trends  occur¬ 
ring  on  the  time  scales  of  days  and  weeks.  The  measured  data  ate  not  a  true 
monthly  average  since  foe  data  ate  rather  sparse  in  time.  Finally,  the  model 
may  not  be  folly  capturing  all  die  oxygen  demand  that  existed  in  this  part  of 
the  Bight,  as  a  result  of  inadequate  depiction  of  loadings,  SOD,  or  some  other 
process. 


In  addition  to  foe  qualitative  model-prototype  comparisons  provided  by 
plots,  a  suite  of  statistical  measurements  were  used  to  compare  observed  and 
simulated  DO  concentrations,  since  DO  is  the  primary  variable  of  interest  in 
this  study.  All  available  DO  measurements  collected  between  15  Apr  and 
30  Sep  1976  were  temporally  and  spatially  associated  with  simulated  predic¬ 
tions  and  the  following  statistics  were  computed: 

1.  Mean  Error,  ME 

ME  *  ~°i)  (1) 

n 

where 

P  =  predicted  value 
O  =  observed  value 
n  =  number  of  observations 
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2.  Absolute  Mean  Enor,  AME 


AME  - 


£(  ft  -  o,  I) 

A 


(2) 


3.  Root-Mean-Square  Enor,  RMS 


RMS  « 


4.  Relative  Error.  RE 

EkI 


(3) 


(4) 


The  calibration  statistics  for  DO  are  given  in  Table  9.  A  total  of  3,664  proto¬ 
type  observations  were  available  for  the  calculation  of  the  statistical  measures. 


Table  9 

CaMbratlon  Statistics  for  DO 

•a** 

Watna 

Maan  Error 

-0.56  gOym4 

Abcoiul*  Maan  Error 

1.69  g  01  m* 

Hoot  Maan-Squaro  Error 

2  20gO,m4 

RMm  Error 

0.09 

A  mean  error  of  aero  is  ideal.  A  negative  mean  error  indicates  that  on  the 
average  the  model  underpredicts,  while  a  positive  mean  error  indicates  over¬ 
prediction.  Examination  of  the  mean  enor  indicates  that  overall  the  model 
underpredicted  measured  values  on  the  average  by  0.55  g  O,  ms.  Although 
some  of  the  discussion  above  focused  on  areas  where  the  model  overpredicted 
DO,  die  overall  underprediction  trend  is  also  evident  in  tire  scatter  (dot  of 
Plate  Dl.  An  absolute  mean  error  of  zero  is  ideal.  The  absolute  mean  error  is 
a  measure  of  the  magnitude  of  the  average  deviations  between  predictions  and 
observations.  The  calculated  absolute  mean  error  of  DO  for  the  calibration 
was  1.69  g  Oj  m  \  The  root-mean-square  error  is  an  indicator  of  the  spread  of 
error  between  predictions  and  observations.  The  relative  error  is  the  ratio  of 
the  absolute  mean  error  to  the  mean  of  the  observations.  The  error  statistics 
are  relatively  high  compared  with  other  modeling  studies,  such  as  Chesapeake 
Bay  (Cerco  and  Cole,  in  preparation),  which  underscores  die  difficulty  of  mod- 


28 


ChapWrS  Appfeaton RssutU 


cling  the  NYB.  The  fact  that  the  sediment  quality  model  was  not  activated 
auy  have  contributed  to  the  error. 


S«n«JtMty  TMts 

Following  preliminary  model  calibration,  the  WQM  was  exercised,  as  dis¬ 
cussed  in  this  section,  to  investigate  Bight  flushing  characteristics  and 
questions  concerning  input  for  model  boundary  conditions.  The  preliminary 
calibration  run  differs  from  the  current  calibration  run  (presented  herein)  in  the 
base  sediment  NH*  release  rate  and  the  temperature  effect  on  SOD.  The  pre¬ 
liminary  calibration  used  a  base  sediment  NH»  rate  of  0.1  g  N  m'2  day*1  and  a 
SOD  temperature  rate  coefficient  of  1.07.  The  preliminary  calibration  run  as 
well  as  the  sensitivity  rest  runs  used  equivalent  kinetic  and  boundary  values 
except  for  the  values  being  rest. 


A  test  was  conducted  to  estimate  "flushing  time”  for  die  NYB.  The  pur¬ 
pose  of  the  test  was  an  order-of-magnitude  evaluation  of  the  time  required  for 
conditions  within  the  Bight  to  be  flushed  by  flows  through  the  boundaries. 

The  flushing  analysis  also  provides  an  estimate  of  die  "memory"  time  of  the 
Bight  The  test  consisted  of  using  a  uniform  initial  conservative  tracer  con¬ 
centration  of  100  units  m*1  in  each  WQM  cell  and  tracer  boundary  condition 
concentrations  of  0.0.  The  boundary  conditions  ensured  that  any  tracer  passing 
out  of  the  Bight  would  not  return.  The  WQM  simulation  extended  from 
IS  Apr  through  30  Sep  1976.  The  number  of  units  remaining  at  the  end  of 
each  month  was  recorded  and  used  to  estimate  the  flushing  time  (Table  10). 
The  number  of  units,  which  is  in  the  form  of  mass,  was  obtained  by  multiply¬ 
ing  the  cell  concentration  times  the  cell  volume. 

Examination  of  Table  10  reveals  that  2.0699  x  10”  units  were  flushed  out 
of  the  NYB  in  154  days  and  that  the  decrease  was  exponential  with  a  decay 
coefficient  of  6.026S  x  10*1  day*1.  Using  the  exponential  decrease  and  the  cal¬ 
culated  decay  coefficient,  50  percent  of  the  tracer  is  flushed  out  of  die  NYB  in 
US  days:  thus,  in  one  year,  only  about  11  percent  of  tire  tracer  is  remaining. 
The  flushing  analysis  therefore  indicates  that  the  flushing  time  of  the  NYB  is 
on  the  order  of  a  year. 

Examination  of  surface  and  bottom  tracer  concentrations  during  the  months 
of  July  and  August  1976  revealed  that  transport  of  tracer  was  upcoast  from  the 
southwest  to  die  northeast  as  evidenced  by  lowest  concentrations  in  the  south¬ 
west  and  highest  concentrations  in  die  northeast  The  simulated  upcoast  circu¬ 
lation  is  consistent  with  observed  hydrographic  flow  anomalies  that  occurred 


TaMiio 

Con— rwttw  Tractr  H-tUng  H— ulf 

Data 

Man  Day 

No.  UaNa  Wawaln—  i  Iff" 

30  Apr 

121 

3.4231 

31  May 

162 

2.5007 

30  Jun 

182 

2.0266 

31  Jui 

213 

1.7827 

31  Aug 

244 

1.5660 

30  Sap 

274 

1.3532 

from  late  spring  until  mid-August  with  the  passage  of  Hurricane  Belle  (NOAA 
1979). 


MocM  boundary  conditions 

Sensitivity  tests  were  conducted  to  examine  the  importance  of  the  benthic 
boundary  (i.e.,  specification  of  SOD)  and  ocean  nutrient  (i.e.,  nitrogen) 
boundary  conditions.  The  interest  in  SOD  stems  from  the  fact  that  the  sedi¬ 
ment  quality  model  was  not  activated  and  SOD  had  to  be  specified  rather  than  •  # 

modeled.  The  primary  data  limitation  exists  along  the  open-ocean  boundaries, 
thus  causing  uncertainty  in  specification  of  ocean  boundary  conditions.  Since 
the  NYB  algae  ar  primarily  nitrogen  limited,  nitrogen  is  a  primary  variable  of 
concern  and  was  the  focus  of  ocean  boundary  sensitivity  tests. 

The  sensitivity  tests  consisted  of  comparing  simulation  results  following  ® 

each  sensitivity  change  with  the  calibration  results.  Comparison  (dots  con¬ 
sisted  of  scatter  [dots,  regional  (e.g..  Apex,  NJ  Nearshore)  vertical  profiles,  and 
transect  plots.  Profiles  and  transects  were  produced  for  June  and  August. 

Scatter  plots  were  for  the  entire  simulation  period  of  15  Apr  through  30  Sep 

1976.  # 

Comparison  statistics  consisted  of  volumetrically  weighted  averages  of  test 
and  calibration  concentrations,  volumetrically  weighted  averages  of  test  and 
calibration  concentration  differences,  and  volumetrically  weighted  RMS  (root 
mean  square)  differences.  Volumetric  weighting  was  employed  to  avoid  hav¬ 
ing  small  and  large  WQM  cells  contribute  equally  to  the  average.  All  compar-  * 

ison  statistics  were  computed  for  the  entire  simulation;  thus,  they  represent  a 
simulation  average.  Monthly  averages  of  concentrations  were  recorded  at  the 
end  of  each  month  for  the  serial  output  of  the  WQM.  The  volume-averaged 
comparison  statistics  were  computed  as  follows: 
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SOD  boundary  condition.  For  the  sensitivity  test  calibration,  a  temporally 
and  spatially  constant  SOD  boundary  condition  of  1.0  g  nr*  da/1  was  used 
because  of  die  limited  spatial  and  temporal  extent  of  the  SOD  data.  A  sensi¬ 
tivity  study  was  conducted  to  evaluate  the  effect  of  inaccurate  SOD  boundary 
conditions.  Two  test  simulations  were  conducted:  (1)  SOD  =  0.0  and  (2)  SOD 
=  10.0  g  m  *  da/1.  The  sensitivity  test  plots  are  displayed  in  Appendix  E  and 
consist  of  the  following: 
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El 


DO  Soft*  PM.  800  -  0.0 


E2 

Vericte  DO  Praflte  Comparison,  Jun,  300  -  00 

ES 

Variate  DO  PioSa  Comparison,  Aug,  800  >  0.0 

E4 

EB 

Vericte  DO  PraMe  Compahaon,  Jun,  800  -  10.0 

ES 

Variate  DO  Praite  Compahaon,  Aug,  800  >  10.0 

The  effect  of  selling  SOD  to  zero  had  a  slight  effect  throughout  the 
modeled  system,  but  more  noticeably  in  the  shallower  areas  of  US,  Apex,  and 
Transect  regions  (Plates  E2  and  E3).  Table  11  lists,  by  region,  simulation- 
average  comparison  statistics  for  the  SOD  =  0.0  test  Examination  of  Table  11 


Tsblsll 

Dissolved  Oxygsn  Comparison  Statistics,  SOD  =  0.0  g  02  m* 
day'1 


1— 

N 

Ctete 

o^or* 

ESEB 

1 

1 

11040 

5.72 

6.27 

0.55 

0.67 

2 

12780 

6.12 

633 

0.12 

0.38 

3 

16120 

5.65 

5.77 

0.12 

0.27 

4 

7380 

5.86 

5.80 

0.04 

008 

6 

7140 

5.48 

5  80 

0.32 

0.44 

6 

5040 

6.00 

6.06 

0.06 

0.14 

7 

4800 

6.25 

6.76 

0.51 

0.63 

8 

22020 

5.80 

6.25 

0.36 

0.50 

9 

10780 

6.06 

6.18 

0.12 

0.23 

10 

4820 

0.54 

10.22 

0.68 

0.02 

11 

10740 

7.63 

8.09 

0.46 

0.53 

NYB 

160060 

6.13 

6.32 

0.10 

0.36 

Note:  8—  TSbls  8  tar  fteyon  nnw.  Region  1 1 1t  US,  and  NYB  i*  Ihe  entire  modeled  sys¬ 
tem.  N  is  tw  number  of  values  used  to  compute  tteiteteSc. 


reveals  that  zero  SOD  sediment  boundary  conditions  increase  the  DO  for  all 
regions  and  by  0.19  g  Oj  m'\  or  3  percent,  for  the  NYB  as  a  whole.  The  most 
notable  change  was  for  Region  10,  the  Transect,  where  DO  increased  7  per¬ 
cent.  Examination  of  Plates  E2  and  E3  (Vertical  DO  Profile  Comparisons) 
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reveals  that  the  most  detectable  increase  in  DO  occurred  in  the  more  shallow 
regions. 


Increasing  the  sediment  boundary  SOD  to  10.0  g  m4  day'1  (an  extremely 
high  and  unrealistic  value)  decreased  the  Bight-wide  DO  by  only  IS  percent 
(Table  12).  However,  this  value  had  a  substantial  effect  in  terms  of  concen¬ 
trations  (Plates  ES  and  E6).  The  most  notable  decrease  was  for  the  Transect 
where  the  DO  decreased  3.S8  g  O,  m  *  (from  9.54  g  O*  m  3  to  5.96  g  Oj  m3), 
a  decrease  of  38  percent  The  high  SOD  also  had  a  sizeable  effect  on 
Regions  7  and  1 1  (LI  Nearshore  and  LIS).  Although  die  value  of  10.0  is 
unrealistically  high,  it  does  demonstrate  some  potential  importance  of  correctly 
specifying  SOD  and  the  need  to  use  the  sediment  model  in  future  work. 


Tabto  12 

Dissolved  Oxygen  Comparison  Statistics,  SOD  =  10.0  g  O,  m  * 
day1 


fr*H. 

N 

Ca m 

UOt** 

Taat 

8  Oi  te’* 

om 

bo^w* 

TOMB  1 

1 

1 

11040 

51.72 

4.00 

-1.63 

1.92 

2 

12780 

6.12 

5.23 

-0.89 

1.31 

3 

15120 

5.65 

5.05 

-0.60 

1.04 

I  4 

7380 

5.85 

5.52 

-0.33 

0.62 

I  5 

7140 

5.48 

4.23 

1.25 

1.62 

1  6 

5040 

6.00 

5.66 

-0.34 

0.70 

I  7 

4800 

6.25 

4.21 

-2.04 

2.33 

1  8 

22920 

5.89 

4.34 

-1.55 

1.95 

9 

18780 

6.06 

5.42 

-0.64 

1.05 

to 

4620 

9.54 

5.96 

-3.58 

4.14 

11 

10740 

7.63 

4.89 

-2.74 

3.02 

NYB 

150060 

6.13 

5.22 

-0.21 

1.44 

Ocean  nitrogen  BC.  A  sensitivity  study  was  conducted  to  investigate  the 
impact  of  variation  in  the  nitrogen  seaward  boundary  conditions.  Two  tests 
were  conducted:  (1)  seaward  boundary  nitrogen  concentrations  set  to  zero 
(N-BC  *  0.0)  and  (2)  an  increase  of  seaward  boundary  nitrogen  concentrations 
by  a  factor  of  two  (N-BC  *  2.0).  Sources  of  seaward  boundary  nitrogen  con¬ 
sisted  of  ammonium-nitrogen,  nitrate-nitrogen,  DON,  PON,  and  algal  nitrogen. 
Algal  nitrogen  sources,  however,  were  not  modified  during  the  sensitivity  tests. 
The  plots  of  jiese  results  are  presented  in  Appendix  E  and  are  identified  as 
follows: 


Chapter  5  AppfeaSon  Raautte 


33 


Plate 

Plot 

E7 

Scatter  Plot,  N-BC  *  0.0 

a 

Vartical  Nanoptankton  Profit*  Comparison,  Jon,  N-BC  *  0.0 

a 

Vortical  Nanoplankton  Profits  Comparison,  Aug,  N-BC  *  0.0 

E10 

Vortical  DO  Prado  Comparison.  Jun,  N-BC  *  0.0 

1  Etl 

Vorttcol  DO  Prado  Comparison,  Aug.  N-BC  *  0.0 

|  E12 

Scatter  Plot,  N-BC  *  2.0 

|  E13 

Vortical  Nanoplankton  Profits  Comparison,  Jun.  N-BC  *  20 

j  E14 

Vortical  Nanoplankton  Profits  Comparison,  Aug,  N-BC  *  2.0 

1  E15 

Vortical  DO  Prods  Comparison,  Jun,  N-BC  *  2.0  | 

i  E16 

Vortical  DO  Prado  Comparison,  Aug.  N-BC  *  2.0  | 

The  time  history  of  monthly  averaged,  global  Bight  nitrogen  concentrations 
is  shown  in  Table  13.  The  flushing  study  indicated  that  the  memory  of  the 
Bight  is  on  the  order  of  a  year  for  a  conservative  substance.  Table  13  indi¬ 
cates  a  shorter  response  time  for  a  nonconservative  substance  such  as  nitrogea 


Table  13 

Bight-Wide  Time  History  for  Total  Nitrogen,  g  N  m* 


Tost 

Ape 

day 

Jun 

Jill 

Aug 

Sap 

|  Calibration 

0.1384 

0.1532 

0.1209 

0.1221 

0.1367 

0.1251 

N-BC  *0.0 

0.0804 

0.0436 

0.0362 

0.0291 

0.0281 

0.0312 

N-BC  *2.0 

0.1067 

0.2612 

0.2066 

0.2162 

0.2469 

0.2209 

The  DO  comparison  statistics  for  the  two  nitrogen  boundary  condition  tests 
are  presented  in  Tables  14  and  IS.  Examination  of  Tables  14  and  IS  reveals 
that  zero  nitrogen  seaward  boundary  conditions  result  in  a  slight  decrease  in 
DO  of  0.S  percent  (Bight-wide)  while  a  doubling  of  the  boundary  conditions 
increases  Bight  DO  by  0.3  percent.  The  WQM  DO  is  relatively  insensitive  to 
the  nitrogen  seaward  boundary  condi tioa  An  increase  in  the  nitrogen  bound¬ 
ary  condition  tends  to  increase  algal  productivity,  which  has  the  effect  of 
slightly  increasing  DO  overall 

Plates  E8  and  E9  (Vertical  Nanoplankton  Profile  Comparisons)  indicate  a 
decrease  in  algae  for  the  lower  nitrogen  boundary  condition.  Nanoplankton 
profiles  are  not  evident  in  Plates  E8  and  E9  for  the  Transect  region  because 
simulated  concentrations  exceeded  the  abscissa  maximum  of  3.0  mg  Chi  m'5. 
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The  corresponding  Vertical  Nanoplankton  Profile  Comparisons  (Plates  E13 
and  E14)  indicate  an  increase  in  algal  mass  with  a  doubling  of  N  boundary 

conditions.  However,  examination  of  Vertical  DO  Profile  Comparisons  with  • 

0.0  nd  doubling  of  the  N  boundary  conditions  (Plates  E10,  El  1,  E1S,  and 

E16)  reveals  minor  deviations  in  die  DO  profiles.  More  impact  cm  DO  could 

be  realized  with  the  sediment  model  activated  since  the  additional  algal  growth 

will  result  in  additional  organic  matter  deposited  on  the  bottom,  thus  increasing 

SOD. 


Demonstration  Scenarios 

The  demonstration  scenarios  include  external  load  increase  and  reduction, 
use  of  die  WQM  for  die  investigation  of  the  cause  of  the  New  Jersey  near¬ 
shore  hypoxia,  and  inclusion  of  LIS  external  loads. 


External  load  Increaaa  and  reduction 

• 

Constant  external  loads  were  varied  for  the  Transect,  NJ  Coast,  and  LI 
Coast  (see  Table  S).  Two  scenario  runs  were  conducted:  (1)  constant  external 
loads  set  to  zero  and  (2)  constant  external  loads  increased  by  a  factor  of  100 
(Table  16).  The  Hudson  River  loads  were  abstracted  from  Table  4.  The  Hud¬ 
son  River,  Sewage  Sludge,  and  Dredged  Material  zone  loadings  were  not 

altered  for  the  scenarios.  Table  17  displays  the  percent  increase  and  decrease  ® 

of  the  loads  for  die  demonstration  scenarios.  Setting  the  constant  external 

loads  to  zero  decreased  the  carbon  and  nitrogen  loads  by  a  factor  of  0.66  and 

0.77,  respectively.  An  increase  of  die  constant  external  loads  by  a  factor  of 

100  increased  die  total  carbon  loads  by  a  factor  of  66  and  the  total  nitrogen 

loads  by  a  factor  of  77.  9 

The  demonstration  scenario  plots  for  external  load  increase  and  decrease  are 
displayed  in  Appendix  F  as  follows: 
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External  Loads,  tons  day* 
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TOC 

nh4-n 

Loads  ■  0 


Tran— ct 

NJ  Coast 

U  Co— t 

Sewape  SM)( 
Oradpad  Malarial 
Hud—n  River 

0. 

0. 

0. 

230. 

1  Constant  External  Loads  *  100 

Transact 

86200. 

NJ  Coast 

7900. 

LI  Co— t 

1200. 

Sowapo  Sludge 
Dredged  Material 
Hud—n  River 

230. 

12600. 

1400. 

5900. 

3900 

600. 

1300. 

300. 

200. 

600 

200. 

200. 
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10. 

7. 

63. 

15. 

52. 

.  19. 
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Table  17 

Percent  Increase  and  Decrease  of  Constant  External  Load 
Scenarios 


POC  NH4-N  |  NOj-N 


External  loads  sat  to  zero 

-81 

-30 

-85 

-36 

-77 

-34 

External  loads  *  100 

7975 

3013 

8382 

3544 

7534 

3879 
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Plate 

PM 

FI 

Scater Plot.  Loads. 0.0 

F2 

Vertical  Nanoplanidon  Pronto  Comparison,  Jon.  Loads  -  0.0 

F3 

Vertical  Nanoptentdon  Profile  Comparison.  Aug,  Loads  -  0.0 

F4 

Vertical  DOC  Profile  Comparison,  Jun,  Loads  •  0.0 

F5 

Vertical  DOC  Profile  Comparison,  Aug,  Loads  »  0.0 

F6 

Vertical  POC  Prefite  Comparison,  Jun.  Loads  «  0.0 

F7 

Vortical  POC  Profile  Comparison,  Aug,  Loads  -  0.0 

F8 

Vortical  DO  Prefite  Comparison,  Jun,  Loads  *  0.0 

F9 

Vertical  DO  Prefile  Comparison,  Aug,  Loads  »  0.0  | 

F10 

Transect  DOC  Comparison,  Jun,  Loads  >  0.0  | 

F11 

Transect  DOC  Comparison,  Aug,  Loads  -  0.0  | 

F12 

Transect  POC  Comparison,  Jun,  Loads  «  0.0 

F13 

Transect  POC  Comparison,  Aug,  Loads  -  0.0 

F14 

Scatter  Plot  Loads  *  100 

F15 

Vertical  Nanoptankton  Profile  Comparison,  Jun,  Loads  *  100 

F16 

Vertical  Nanopianklon  Profile  Comparison,  Aug,  Loads  *  100 

F17 

Vertical  DOC  Profile  Comparison,  Jun,  Loads  *  100 

F18 

Vertical  DOC  Profile  Comparison,  Aug,  Loads  *  100 

F18 

Vertical  POC  Profile  Comparison,  Jun,  Loads  *  100 

F20 

Vertical  POC  Profile  Comparison,  Aug,  loads  *  100 

F21 

Vertical  DO  Profile  Comparison,  Jun,  Loads  *  100 

F22 

Vertical  DO  Profile  Comparison,  Aug,  Loads  *  100  | 

F23 

Transect  OOC  Comparison,  Jun,  Loads  *  100  | 

F24 

Transect  DOC  Comparison,  Aug,  Loads  *  100 

F25 

Transect  POC  Comparison,  Jun,  Loads  *  100 

F26 

1 

Transect  POC  Comparison,  Aug,  Loads  *  100 

In  addition,  comparison  statistics  for  simulation-averaged  DO  were  computed 
regional  and  Bight-wide  (i.e.,  NYB)  and  are  provided  in  Tables  18  and  19  for 
loads  =  0.0  and  loads  *  100,  respectively. 

Decreasing  the  external  loads  had  the  effect  of  decreasing  algal,  DOC,  and 
POC  concentrations.  Although  organic  carbon  decreased,  the  net  effect  of 
these  changes  resulted  in  a  slight  decrease  in  DO,  most  noticeable  in  the  Tran¬ 
sect  region.  The  decrease  in  algae  had  more  impact  on  decreasing  DO  than 
die  effect  of  decreasing  organic  carbon  had  on  increasing  DO.  Dramatic 
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loading  increases  caused  slight  DO  decrease  Bight-wide  and  a  substantial  DO 
decrease  in  the  Transect  region.  Nearshore  regions  (e.g.,  die  Apex)  also  exhib¬ 
ited  substantial  decreases.  The  loading  increase  also  caused  high  algal  and 
organic  carbon  concentrations. 

Overall,  the  whole  NYB  system  exhibits  little  sensitivity  to  loading 
changes,  especially  for  IX)  (see  NYB  in  Tables  18  and  19),  demonstrating  that 
the  NYB  has  a  tremendous  assimilative  capacity.  The  above  demonstrations 
should  not  be  interpreted  as  what  would  happen  with  loading  changes.  It  is 
pointed  out  that  this  model  does  not  include  the  sediment  model  that  can 
accumulate  organic  matter  and  adjust  the  SOD  accordingly.  These  simulations 
merely  demonstrate  how  the  model  might  be  used  and  roughly  estimate  the 
sensitivity  of  the  Bight  to  loading  changes.  For  more  definitive  estimates  of 
the  effects  of  loading  changes  on  water  quality,  the  sediment  model  should  be 
activated  and  calibrated,  and  long-term  (e.g.,  multiyear)  simulations  would  be 
required  to  properly  capture  sedimentation  and  mineralization  of  organic  matter 
in  the  sediments  and  the  resulting  SOD.  Additionally,  the  WQM  should  be 
more  accurately  calibrated  for  longer  periods  of  time. 

The  external  loads  were  inserted  in  cells  immediately  adjacent  to  shore. 
However,  both  the  load  comparison  statistics  and  die  plots  used  spatially 
averaged  results  (averaged  over  regions  and  the  Bight).  Near-field  effects  that 
may  occur  within  meters  to  a  few  kilometers  of  the  point  of  discharge  would 
be  attenuated  through  the  averaging  process  and  may  not  be  evident  in  the 
regional  averages. 


1  Table  18 

|  Dissolved  Oxygen  Comparison  Statistics,  Loads  =  0.0 

Region 

N 

CaNb 

gOjie'* 

Teat 

g^nf* 

DM 

g  Oj  in'* 

TRMS 

80**-* 

1 

11040 

5.72 

5.66 

-0.06 

0.13 

2 

12780 

6.12 

6.12 

-0.00 

001 

3 

15120 

5.65 

5.65 

-0.00 

0.00 

4 

7380 

5.85 

5.85 

0.00 

0.00 

S 

7140 

5.48 

5.48 

-0.00 

0.00 

6 

5040 

6.00 

6.00 

-0.00 

0.00 

7 

4800 

6.25 

6.23 

-0.02 

0.04 

8 

22820 

588 

5.88 

-0.00 

0.02 

8 

18780 

6.06 

6.06 

-0.00 

0.00 

10 

4820 

8.54 

7.70 

-1.84 

2.82 

11 

10740 

7.63 

7.54 

-0.06 

0.24 

NYB 

150060 

6.13 

6.12 

-0.00 

0.08 
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TaMa  19 

ntaaoftratl 

OxyganCc 

Mnpanmn  % 

Lateda  *  100 

P— 

N 

•<^«4 

ns 

^391 

TMM 

9  ■"* 

1 

11040 

5.72 

5.14 

-0.58 

1.33 

2 

12780 

8.12 

6.13 

0.01 

0.16 

a 

18120 

5.65 

5.66 

•0.00 

0.06 

4 

7380 

5.86 

6.84 

-0.00 

0.02 

6 

7140 

6.48 

5.42 

-0.06 

0.37 

6 

5940 

6.00 

5.96 

-0.01 

0.12 

7 

4800 

6.25 

5.50 

-0.66 

1.14 

8 

22920 

5.89 

5.66 

-0.23 

0.77 

9 

18780 

6.06 

6.02 

-0.03 

0.37 

10 

4820 

9.54 

2.20 

-7.34 

7.96 

11 

10740 

7.63 

7.67 

0.04 

0.31 

NYB 

150060 

6.13 

6.03 

-0.10 

0.52 

New  Jersey  neershore  hypoxia  • 

The  measured  bottom  layer  DO  contour  shown  in  Plate  23,  which  was 
obtained  from  the  bottom  5-m  oxygen  distribution  presented  in  Stoddard 
(1983,  1989),  indicates  that  low  DO  occurred  during  July  1976  off  the  coast  of 

New  Jersey.  Stoddard’s  model  oxygen  budget  investigations  indicated  that  . 

transport  was  a  key  factor  in  the  onset  of  hypoxia.  The  onset  of  steady  south¬ 
west  winds  in  latter  May  resulted  in  an  upcoast  northeast  flow  with  upwelling. 

The  upcoast,  upwelling  flows  triggered  the  accumulation  of  Ceratium  off  the 
coast  of  New  Jersey  whose  respiration  and  decay  contributed  to  oxygen  deple¬ 
tion.  However,  the  NYB  WQM  simulation  predicted  small  biomass  of  net 

plankton,  which  was  not  in  disagreement  with  small  measured  values.  • 

Although  few  measured  data  were  available  for  net  plankton,  their  magnitude 
suggested  small  concentrations  for  initial  and  ocean  boundary  conditions. 

The  WQM  was  used  to  investigate  why  a  DO  minimum  was  correctly 
simulated  (see  Plate  23,  where  "Simulated”  results  were  obtained  from  the  m 

WQM  of  this  present  study)  off  the  coast  of  New  Jersey.  Plate  26  displays  the 
4.0-g  O2  ni  3  contour  plot  for  the  calibration  ran.  The  4.0-g  02  m'3  contour 
was  selected  for  the  New  Jersey  nearshore  hypoxia  investigation  because  this 
contour  appeared  in  all  the  test  investigations. 

The  first  step  in  the  investigation  was  to  examine  the  magnitude  of  foe  • 

source  and  sink  terms  in  the  DO  kinetic  terms.  The  DO  kinetic  terms  were 
recorded  for  foe  water  column  corresponding  to  WQM  cell  1123  (ij  =  16,17). 
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WQM  cell  1125  was  selected  because  it  occurred  in  the  center  of  the  July 
minimum.  The  water  column  depth  was  31.87  m.  The  kinetic  DO  fluxes  for 
each  layer  were  typically  negative.  The  DO  kinetic  fluxes  indicated  that  DOC 
mineralization  was  the  dominant  kinetic  DO  sink  for  layers  1-9  and  SOD  was 
the  dominant  sink  for  layer  10,  the  bottom  layer. 

An  examination  of  the  kinetic  sources  and  sinks  of  DOC  indicated  that  foe 
mineralization  sink  of  DOC  greatly  exceeded  algal  sources  and  sources  through 
hydrolysis  of  POC.  Alternative  sources  of  DOC  include  advection  of  DOC 
from  foe  New  Jersey  coastline  (NJ  external  constituent  loads  listed  in  Table  5) 
and  advection  from  the  seaward  boundary.  The  WQM  does  provide  options 
for  recording  advective  fluxes  through  cell  walls,  but  identification  of  the 
source  required  examination  of  the  NJ  coast  external  loads  and  boundary 
conditions. 

A  WQM  simulation  with  NJ  coast  external  loads  set  to  zero  resulted  in  no 
observed  difference  in  DO  or  DOC  in  the  water  column  corresponding  to 
WQM  cell  1 125  or  the  spatial  pattern  of  foe  4.0-g  02  m  3  contour  (Plate  27). 
The  northeast  flows  along  foe  New  Jersey  coast  did  not  advect  the  coastal 
loads  to  WQM  cell  1 125.  The  remaining  source  of  DOC  was  the  ocean 
boundary  conditions.  The  northeast  flows  suggest  that  the  southwest  boundary 
conditions  for  Segments  1,  2,  and  3  (Plate  5)  could  possibly  affect  WQM  cell 
1 125.  A  WQM  simulation  was  conducted  with  Segments  1,  2,  and  3  DOC 
boundary  concentrations  set  to  0.0.  A  comparison  of  DO  contour  {dots  for  the 
calibration  run  (Plate  26)  and  DOC  boundary  conditions  for  Segments  1,  2, 
and  3  set  equal  to  zero  (Plate  28)  revealed  that  foe  patterns  of  tow  DO  concen¬ 
trations  were  similar.  However,  the  spatial  extent  of  low  DO  was  smaller  with 
the  DOC  boundary  condition  segments  set  to  zero.  The  flux  of  DOC  from  the 
ocean  boundaries  intensified  DO  depletion  off  the  coast  of  New  Jersey,  but  it 
was  not  foe  major  cause  of  foe  hypoxia. 

Because  SOD  was  the  dominant  kinetic  sink  for  DO  for  the  bottom  layer, 
a  WQM  simulation  was  conducted  with  SOD  set  equal  to  0.0  for  the  NJ  Near¬ 
shore,  NJ  Midshelf,  rod  NJ  Shelf  regions.  The  pattern  of  the  simulated  DO 
contour  off  the  New  Jersey  coast  is  substantially  different  with  the  SOD  set  to 
zero  (Plate  29),  indicating  that  SOD  played  an  important  role  in  the  New  Jer¬ 
sey  nearshore  hypoxia.  However,  there  is  still  a  relatively  large  zone  of  low 
DO  with  the  SOD  set  to  zero,  indicating  other  remaining  causes  of  foe 
hypoxia. 

The  DO  ocean  boundary  conditions  were  examined.  The  DO  boundary 
concentrations  for  Segments  1,  2,  and  3  listed  in  Appendix  C  are  reproduced  in 
Table  20. 

The  significance  of  the  DO  boundary  condition  was  investigated  by  reflect¬ 
ing  foe  DO  boundary  condition.  Reflection  of  the  DO  boundary  condition 
means  that  foe  boundary  DO  concentration  was  set  equal  to  the  immediately 
interior  WQM  cell  DO  concentration  from  foe  previous  time-step.  The  DO 
boundary  conditions  were  thus  completely  dependent  upon  transport  and 
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kinetics  interior  to  the  boundary.  Examination  of  contour  (dots  of  bottom  layer 
DO  for  July  with  reflective  DO  boundaries  revealed  that  DO  was  greater  than 
4.0  g  Oj  m'3  globally  (i.e..  Bight-wide).  In  contrast,  the  calibration  bottom 
DO  during  July  was  generally  less  that  4.0  off  the  coast  of  New  Jersey  with 
extensive  areas  of  DO  less  than  3.0.  The  conclusion  of  the  reflective  DO 
boundary  experiment  was  that  low  DO  concentrations  for  the  southwest  ocean 
boundary  were  a  major  contributor  to  simulated  DO  depression  off  the  coast  of 
New  Jersey. 

The  New  Jersey  nearshore  hypoxia  investigation  revealed  that  low  DO 
simulated  off  the  coast  of  New  Jersey  was  die  result  of  three  major  interacting 
components:  (1)  the  prevailing  southwest  to  northeast  residual  flows,  (2)  DOC 
and  DO  boundary  conditions  along  the  southwest  ocean  boundary  of  die  model 
grid,  and  (3)  SOD.  The  prevailing  southwest  to  northeast  residual  flows 
advected  the  southwest  boundary  condition  constituents  upcoast  along  the  New 
Jersey  coast  The  southwest  boundary  provided  a  source  of  DOC  of  larger 
magnitude  than  from  algae  and  hydrolysis  of  POC.  The  southwest  boundary 
condition  for  DO  was  specified  using  measured  prototype  data.  Measured 
bottom  DO  used  in  the  specification  of  the  southwest  boundary  conditions  was 
frequently  less  than  3.0  g  0^  m  3.  The  prevailing  residual  currents  advected 
the  low  DO  boundary  waters  along  the  NJ  coastline.  The  mineralization  of  the 
DOC  intensified  DO  depletion.  The  SOD  increased  die  rate  of  DO  depletion 
in  the  bottom  layer.  This  increased  rate  of  DO  depletion  was  manifested 
through  earlier  occurrence  and  greater  intensity  of  low  DO  concentrations. 

The  New  Jersey  Coast  external  constituent  loads  had  little  influence  on  the 
simulated  DO  depression  off  the  coast  of  New  Jersey. 


Chapter  5  Application  Fteautte 


As  staled  earlier,  calibrations,  sensitivity  tests,  and  load  reduction  scenarios 
were  conducted  without  LIS  external  loads.  An  investigation  was  conducted  to 
examine  the  impact  of  US  external  loads  on  the  DO  resources  of  the  Bight 

The  US  external  loads  woe  based  on  tabulations  listed  in  Wolfe  et  al. 
(1991)  (Table  21). 


Tab*  21 

Coolant  External  Conatttuant 1  pads,  LIS,  tons  day'1 

DOC 

roc 

MVN 

HOfH 

DON 

FOti 

411. 

137. 

101. 

16. 

10. 

Table  22  lists  by  region  die  simulation-average  DO  comparison  statistics  for 
the  US  external  load  test 


Table  22 
Diaaolvad 
External  I 

Oxygen  Cc 
joada 

Mnpartaon  Statistics,  L 

ang  Island  Sound  j 

BsAin 

nayaNi 

N 

Cafe 

ap*".-* 

T«at 
a  Of"1"* 

DM 

man 

ac^w-* 

1  i 

11040 

606 

6  06 

0.00 

0.01 

2 

12780 

6.30 

6.30 

0.00 

000 

3 

15120 

5.87 

5.87 

0.00 

nrm 

0.00 

4 

7380 

6.02 

6.02 

0.00 

0.00 

5 

7140 

5.88 

5.88 

0.00 

0.00 

6 

5040 

6.13 

6.13 

0.00 

0.00  | 

7 

4800 

669 

6.69 

0.00 

0.00  1 

8 

22B20 

6.33 

6.33 

0.00 

000  i 

0 

18780 

6.28 

6.29 

0.00 

0.00 

10 

4620 

9.71 

9.82 

0.10 

0.16 

11 

10740 

7.87 

7.99 

0.11 

0.20 

NYB 

150080 

6.40 

6.40 

0.00 

0.03 

Examination  of  Table  22  reveals  that  inclusion  of  LIS  external  loads 
increases  the  regional  average  DO  of  the  Transect  region  0.10  g  02  m'3  and  of 
the  US  region  0.11  g  02  m  .  A  difference  in  DO  was  not  detected  in  other 
regions  or  Bight-wide.  Examination  of  the  DO  vertical  profile  comparison  for 
August  (Plate  30)  indicates  the  small  increase  in  DO  in  die  Transect  and  US. 
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An  additional  conclusion  from  the  examination  of  Plate  30  is  that  the  simu¬ 
lated  DO  in  US  was  uniformly  distributed  from  surface  to  bottom.  The  uni¬ 
form  distribution  may  be  due  to  excessive  vertical  diffusion  and/or  inability  to 
capture  vertical  stratification  in  US  using  12.5-hr  intertidal  hydrodynamic 
averaging. 


The  slight  increase  in  DO  was  the  result  of  an  increase  in  algal  photosyn¬ 
thesis,  a  pattern  similar  to  the  sensitivity  response  due  to  increased  nitrogen 
boundary  concentrations.  Dissolved  oxygen  depletion  associated  with 
increased  algal  growth  followed  by  mortality,  settling,  and  decay  of  algal  car¬ 
bon  in  the  sediments  is  not  simulated  without  the  inclusion  of  the  sediment 
diagenesis  model.  The  DO  hypoxia  in  die  Western  Basin  of  die  Sound  (Parker 
and  O’Reilly  1991)  was  not  captured  by  the  model.  Possible  reasons  that  the 
hypoxia  was  not  simulated  are  (1)  the  loads  were  uniformly  distributed  from 
East  River  to  Montauk  Point  and  (2)  excessive  vertical  mixing  prevented  phys¬ 
ical  stratification. 
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Recommendations 


Conclusions 

This  study  demonstrates  that  a  water  quality  model  of  the  NYB  is  feasible 
and  provides  reasonable  results.  The  model  simulates  the  Bight,  LIS,  harbor, 
and  estuaries  as  a  single,  dynamically  coupled  unit  This  approach  has  much 
merit  for  assessing  how  the  various  regions  interact  and  for  providing  a  true 
representation  of  die  real  system. 

A  model  such  as  this  provides  a  strong  capability  for  examining,  in  a  cost- 
effective  manner,  a  host  of  questions  that  could  not  otherwise  be  addressed. 
The  present  model  is  not  developed  to  the  point  that  it  can  be  used  for  deter¬ 
mining  future  nutrient  loading  goals  or  waste  load  allocations  (WLAs),  but  it  is 
a  large  step  toward  such  a  model,  which  is  now  considered  attainable,  hi  its 
present  state  of  development,  the  model  is  well  suited  to  examining  a  number 
of  questions,  some  of  which  are  discussed  herein. 

The  model  compared  relatively  well  with  observations  in  die  Bight  and 
successfully  captured  die  summer  hypoxia  of  1976.  However,  witter  quality 
conditions  within  the  Transect  region  (i.e.,  in  the  harbor  and  estuaries)  did  not 
compare  well,  and  improvements  are  needed  if  any  additional  modeling  is 
undertaken.  Recommendations  for  these  improvements  are  discussed  below. 
The  primary  reasons  for  the  hypoxia  off  the  New  Jersey  coast  during  die  sum¬ 
mer  of  1976,  as  determined  by  die  model,  were  SOD  and  the  advectkm  of 
DOC  and  low  DO  from  the  southwest  ocean  boundary  toward  the  northeast. 
The  advecdon  of  low  DO  had  the  greatest  effect  of  these  three  processes. 

Loading  sensitivity  studies  indicated  that  throughout  most  of  the  Bight, 
water  quality  (e.g.,  algae  and  DO)  is  insensitive  to  loading  decreases,  but  the 
Transect  region  and  part  of  the  Apex  are  sensitive.  Extremely  large  loading 
increases  were  required  to  register  effects  in  the  Bight  Although  algae,  POC, 
and  DOC  increased  substantially  (especially  closer  to  shore  and  within  the 
Transect)  when  total  loading  was  increased  over  60-fold,  DO  within  die  Bight 
changed  little.  The  effect  on  SOD  of  this  loading  increase  was  not  accounted 
for  since  the  sediment  diagenesis  model  was  not  activated.  However, 
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increasing  the  SOD  to  10  g  mJ  day'1  (an  unrealistically  high  value)  only  mod¬ 
erately  reduced  DO  in  Bight  bottom  waters  (e.g.,  a  reduction  of  about  1  to  3  g 
DO  m1).  It  appears  that  the  tremendous  volume  of  the  Bight  provides  a  great 
assimilative  opacity  for  DO,  which  tends  to  increase  in  the  seaward  direction. 
Setting  the  SOD  to  zero  only  slightly  increased  bottom  water  DO  by  about 
1.0  g  DO  m'1  in  nearshore  regions,  and  even  less  farther  offshore.  Sensitivity 
tests  on  the  ocean  boundary  condition  for  nitrogen  indicated  that,  although 
algae  responded,  DO  of  the  NYB  system  was  insensitive  to  changes  in  nutri¬ 
ents  introduced  from  the  ocean. 


Modeling  Recommendations 

The  model  results  point  out  several  deficiencies  that  should  be  overcome  if 
the  model  is  to  be  used  for  establishing  nutrient  reduction  goals  and  WLAs. 
Several  steps  should  be  taken  to  obtain  a  more  accurate  WQM  as  discussed 
below.  These  recommendations  are  not  necessarily  in  order  of  priority. 

a.  The  sediment  diagenesis/flux  model  should  be  activated  and  calibrated. 
Recent  sediment  flux  measurements  from  LIS  could  be  used  to  assist  in 
this  calibration.  Activation  of  the  sediment  model  may  also  require 
activation  of  additional  water  quality  state  variables. 

b.  Additional  data  sources  should  be  sought  for  better  defining  loadings 
and  boundary  conditions  and  for  making  calibration  comparisons.  Insuf¬ 
ficient  water  quality  data  were  utilized  in  this  study  for  foe  Transect 
region  and  LIS,  which  makes  it  difficult  to  draw  meaningful  conclusions 
about  model  accuracy  in  these  areas.  More  accurate  loading  estimates 
for  the  Hudson  River  should  also  be  sought 

c.  More  frequent  update  intervals  for  loadings  and  boundary  conditions 
should  be  used.  For  example,  daily  or  weekly  Hudson  River  flows  and 
concentrations  (or  loads)  should  be  used.  At  least  monthly  updates 
should  be  used  for  external  loads. 

d.  The  excessive  algal  concentration  computed  within  the  Transect  region 
should  be  investigated  further  and  remedied.  The  most  probable  solu¬ 
tion  to  this  problem  is  to  include  the  effect  of  suspended  solids  mi  light 
attenuation.  If  there  are  insufficient  data  to  include  suspended  solids  as 
a  model  state  variable,  it  may  be  possible  to  indirectly  relate  light  atten¬ 
uation  to  suspended  solids  through  a  relationship  to  Hudson  River  flow 
discharge. 

e.  Refractory  DOC  can  be  added  to  the  model  to  better  match  measured 
DOC. 
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/.  Multiyear  (e.g.,  2  years  or  more)  calibntkw/verifkadon  simulations 
should  be  conducted  to  obtain  a  more  accurate  model  and  more  fully 
demonstrate  the  accuracy. 

g.  A  primary  deficiency  concerns  the  difficulty  in  representing  the  proper 
vertical  mixing.  During  the  course  of  the  study,  the  parameters  that 
affect  the  vertical  eddy  diffusivity  calculations  in  the  HM  were  varied, 
causing  WQM  results  for  temperature  and  DO  to  range  from  well-mixed 
to  overiy  stratified.  A  higher  order  closure  algorithm  for  vertical  eddy 
diffusivity,  such  as  transport  equations  for  kinetic  energy  and  dissipa¬ 
tion,  may  be  required  to  properly  simulate  vertical  mixing. 

There  are  at  least  two  investigations  that  should  be  undertaken  to  improve 
confidence  in  the  model.  These  investigations  may  lead  to  additional  model 
improvements  or  modifications  as  explained  below. 

a.  Investigate  transport  in  more  detail.  Vertical  diffusion  is  computed  as  a 
simple  arithmetic  mean  in  the  intertidal  averaging  procedure.  This  pro¬ 
cedure  may  not  be  adequate  in  the  more  dynamic  regions  of  the  harbors, 
estuaries,  and  LIS.  The  proper  preservation  of  vertical  mixing  with 
intertidal  transport  was  not  investigated  in  this  study.  This  was  found 
not  to  be  a  problem  in  the  Chesapeake  Bay  model  study  (Dortch  1990), 
but  that  estuary  is  weakly  nonlinear  whereas  the  inner  regions  of  this 
study  site  are  much  more  nonlinear.  In  any  future  modeling  of  this 
system,  this  aspect  should  be  studied  more  closely.  It  may  be  possible 
to  still  use  intertidal  currents,  but  the  vertical  diffusion  may  have  to  be 
handled  differently.  Additionally,  tests  should  be  conducted  to  ensure 
that  mass  flux  through  the  Transect  is  property  represented  with  inter¬ 
tidal  transport.  This  test  can  be  accomplished  by  comparing  HM  mass 
or  salt  flux  through  the  Transect  with  that  of  the  WQM  using  intertidal 
hydrodynamic  updates. 

b.  Conduct  additional  investigations  of  the  sensitivity  of  ocean  boundary 
conditions.  Further  development  of  the  model  for  determining  nutrient 
reduction  goals  would  subsequently  involve  simulating  various  proposed 
nutrient  reduction  strategies  to  evaluate  their  effectiveness.  As  dismissed 
earlier,  the  investigation  of  the  cause  of  the  hypoxia  revealed  the  impor¬ 
tance  of  ocean  boundary  conditions  for  DO  and  DOC.  This  brings  up  a 
question:  How  can  water  quality  with  future  loading  conditions  be  pre¬ 
dicted  if  the  results  depend  on  ocean  boundary  conditions  which  could 
depend  on  the  future  loadings?  Fortunately,  die  results  of  this  study 
generally  support  the  idea  that  loadings  from  the  land  boundaries  have 
minimal  effect  on  conditions  near  the  ocean  boundaries.  However,  this 
idea  should  be  more  thoroughly  investigated  to  ensure  that  nutrient 
loading  strategies  are  simulated  property. 

If  future  modeling  must  address  WLAs  and  their  impacts  within  die  Tran¬ 
sect  region,  a  fine-grid  model  of  the  inner  harbors  and  bays  is  recommended. 

It  is  not  possible  to  accurately  resolve  conditions  inside  the  Transect  with  the 
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present  grid.  The  simplest  approach  would  be  to  construct  a  separate  model  of 
the  inner  harbors  and  bays  which  extends  into  the  Apex  and  LIS  so  that  it  can 
be  interfaced  with  the  NYB  model  (i.e.,  the  grid  of  the  present  model).  The 
models  would  not  have  to  be  dynamically  coupled;  rather,  results  from  the 
NYB  model  could  be  used  to  set  boundary  conditions  for  the  fine-scale 
haibor/estuaries  model. 


Water  Quality  Monitoring  Recommendations 

This  section  provides  general  recommendations  for  water  quality  monitoring 
for  two  purposes;  to  support  a  future  comprehensive  eutrophication  model  and 
to  monitor  future  water  quality/eutrophication  trends.  The  first  purpose 
requires  detailed,  synoptic  sampling  over  a  rather  short  time  frame,  while  the 
second  purpose  requires  much  less  detailed  sampling  over  a  long  time  frame. 
These  two  approaches  to  monitoring  are  referred  to  here  as  synoptic  monitor¬ 
ing  for  modeling  and  long-term  monitoring.  Both  types  of  monitoring  are  for 
general  water  quality  (e.g.,  nutrients  and  DO)  and  do  not  address  contaminants 
(e.g.,  trace  metals,  hydrocarbons,  and  synthetic  organics). 

Neither  monitoring  purpose  would  be  used  for  real-time  water  quality  mod¬ 
eling.  Real-time  modeling  is  defined  here  as  trying  to  predict  conditions  in  the 
near  future  based  on  recent  observations.  Eutrophication/general  water  quality 
models  of  large  water  bodies  such  as  this  are  not  used  for  real-time  water 
quality  prediction,  but  rather  for  planning  and  evaluating  future  management 
strategies  (e.g.,  evaluating  the  effects  of  reducing  nutrient  loadings).  The  water 
quality  observed  at  a  particular  location  may  be  the  result  of  months  (even 
years)  of  previous  conditions  because  of  the  long  memory  of  large  water  bod¬ 
ies.  For  example,  low  DO  observed  during  August  may  be  the  result  of  load¬ 
ings  and  algal  blooms  that  occurred  in  the  spring.  Therefore,  it  is 
unreasonable  to  expect  such  water  quality  models  to  predict  what  the  water 
quality  will  be  next  month  based  on  the  observed  water  quality  this  month. 
With  this  in  mind,  water  quality  monitoring  should  not  be  conducted  to  sup¬ 
port  real-time  water  quality  modeling. 

The  discussion  below  provides  general  guidance  for  conducting  the  two 
types  of  monitoring  programs.  The  details  for  these  monitoring  programs  are 
not  provided  since  formulation  of  such  details  would  depend  on  logistical 
capabilities  of  the  parties  that  would  collect  these  data  and  other  conditions 
discussed  below.  Additionally,  the  effort  required  to  formulate  the  details  of 
the  monitoring  programs  is  beyond  the  scope  of  this  study. 


Synoptic  monitoring  for  modeling 

The  water  quality  model  feasibility  study  described  herein  was  constrained 
because  of  data  limitations.  The  constraint  was  in  part  based  on  the  study 
period  selected,  i.e.,  DO  hypoxia  event  during  summer  1976.  There  may  be 
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yean  where  more  complete  data  were  obtained  (such  as  monitoring  of  addi¬ 
tional  nutrient  forms).  If  a  comprehensive  eutrophication  model  is  performed, 
more  complete  data  will  be  required  to  ensure  proper  model  implementation. 
Such  data  sets  are  referred  to  as  synoptic  and  would  be  used  primarily  for 
water  quality  model  calibration  and  verification. 

Before  a  synoptic  water  quality  monitoring  program  is  developed  and  exe¬ 
cuted.  the  existing  database  should  be  examined  more  closely.  Stoddard  and 
Han  (1989)  reviewed  rod  assessed  field  data  available  for  use  in  this  study. 
Their  inventory  of  available  data  indicated  that  an  enormous  amount  of  water 
quality  data  was  readily  available  from  databases  to  characterize  the  water 
quality  of  die  New  York  Bight.  They  provided  WES  with  data  from  the 
period  of  study  (i.e.,  1976).  Relatively  complete  data  sets  may  exist  for  other 
periods  as  well,  in  addition  to  more  recently  collected  data,  e.g.,  sediment  flux 
data  collected  in  Long  Island  Sound. 

The  first  step  in  a  comprehensive  eutrophication  model  study  should  be  to 
collate  and  study  the  existing  database  and  to  select  study  years  where  the  data 
are  most  complete.  Data  gaps  should  be  identified  to  help  determine  what,  if 
any,  additional  data  are  required  for  a  comprehensive  eutrophication  model 
study.  Problems  may  be  encountered  with  the  historical  data,  such  as  changes 
in  collection  and  analytical  techniques  that  occurred  over  foe  last  two  decades. 
If  the  existing  data  are  deemed  adequate,  a  synoptic  monitoring  program  is  not 
needed.  However,  if  the  existing  data  are  considered  inadequate  to  support 
such  a  modeling  effort,  a  synoptic  data  set  based  on  current  Bight  conditions 
and  analytical  techniques  should  be  considered  as  described  below. 

A  3-year  data  collection  program  would  be  adequate  for  establishment  of  a 
synoptic  data  set  The  constituents  measured  should  include  salinity,  tempera¬ 
ture,  chlorophyll  and  algal  quantification  and  speciation,  multiple  forms  of 
organic  carbon,  multiple  forms  of  organic  and  inorganic  nitrogen  and  phospho¬ 
rus.  multiple  forms  of  silica,  and  dissolved  oxygen.  Additionally,  some  work 
on  partitioning  of  labile  and  refractory  organic  components  is  recommended. 

Year-round,  bimonthly  sampling  stations  should  be  located  in  the  Hudson, 
Raritan,  and  Passaic  Rivers  just  upstream  of  their  fall  lines.  Bimonthly 
sampling  should  also  be  conducted  between  April  and  October  for  stations  in 
Long  Island  Sound,  Newark  and  Raritan  Bays,  foe  Kills,  and  foe  Sandy  Hook- 
Rockaway  transect.  Monthly  sampling  should  be  conducted  at  these  stations 
during  winter  months.  Three  to  six  stations  are  needed  in  Long  Island  Sound. 
Newark  and  Raritan  Bays  and  the  Kills  should  have  one  or  two  stations  each. 
Approximately  three  stations  should  be  located  across  foe  transect 

Monthly  sampling  should  be  conducted  on  the  seaward  boundaries  with  a 
station  for  each  of  the  seven  boundary  segments  shown  in  Plate  21.  Monthly- 
to-seasonal  (shout  every  2  months)  sampling  should  be  conducted  in  foe 
remaining  extent  of  the  Bight  at  approximately  nine  stations  which  represent 
the  nine  regions  outside  foe  transect  in  Plate  23.  For  these  nine  Bight  stations, 
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monthly  sampling  should  occur  hum  April  through  October  with  sampling 
every  2  months  during  the  winter. 

All  sampling,  except  the  fail  lines,  should  be  conducted  over  the  water  col¬ 
umn,  i.e.,  at  multiple  depths  for  each  station.  A  minimum  of  two  depths  (i jt., 
near-surface  and  near-bottom)  are  required,  but  additional  depths  between  these 
two  are  highly  desirable,  such  as  within  the  pycnochne. 

A  comprehensive  eutrophication  model  study  will  require  activation  of  the 
sediment  di  agenesis/flux  model  The  sediment  model  oxygen  and  nitrogen 
parameters  measured  for  Chesapeake  Bay  were  transferable  to  Long  Island 
Sound  and  Namgansett  Bay  applications;  however,  the  phosphorus  parameters 
were  not  transferable.  Sediment  water  fluxes  of  oxygen,  nitrogen,  phosphorus, 
silicon,  and  sulfide  should  be  measured  for  the  Bight  at  multiple  locations.  A 
sediment  sampling  station  should  be  established  in  each  of  the  major  interior 
areas,  i.e.,  Hudson  River,  Raritan  Bay,  Newark  Bay,  the  Kills,  and  within  the 
transect  region  (see  Plate  23).  Nine  sites  should  also  be  located  throughout  the 
Bight  to  cover  the  nine  regions  shown  outside  the  transect  in  Plate  23.  Sedi¬ 
ment  sampling  for  flux  measurements  should  be  conducted  four  times  a  year 
(winter,  spring,  summer,  and  fall)  throughout  the  3-year  synoptic  monitoring 
study.  Each  sample  should  be  extracted  as  an  intact  core  of  about  30-cm 
length.  The  flux  measurements  should  be  conducted  immediately  onboard  the 
survey  ship  following  collection.  The  same  procedures  used  in  foe  Chesapeake 
Bay  mid  Long  Island  Sound  studies  should  be  followed. 

The  above  synoptic  monitoring  program  would  cost  several  million  dollars. 
This  is  not  a  firm  cost  proposal,  but  rather  a  rough  estimate  for  planning  pur¬ 
poses.  This  work  would  most  likely  be  contracted  to  an  organization  that  does 
estuarine  and  oceanographic  monitoring  routinely  and  has  foe  capability  to 
cover  extensive  geographic  areas  in  a  relatively  short  time  frame. 

It  is  emphasized  that  the  final  monitoring  recommendations  would  depend 
on  foe  results  of  a  data  compilation  and  synthesis  effort  and  budgetary  con¬ 
straints.  The  final  recommendations  could  extend  from  initiating  few,  if  any, 
of  the  items  in  foe  above  plan  to  all  of  foe  items  in  foe  above  plan. 


Long-term  monitoring 

A  separate  and  very  different  monitoring  {dan  would  be  used  to  trade  long¬ 
term  water  quality  trends.  Water  quality  is  also  constrained  here  to  nutrients, 
DO,  and  phytoplankton  (i.e.,  eutrophication).  This  type  of  data  would  be  used 
to  monitor  foe  health  of  foe  Bight  and  how  that  health  might  be  changing  with 
time.  Such  data  would  not  be  fully  adequate  for  modeling  future  water 
quality/eutrophication. 

Several  monitoring  programs  are  ongoing,  as  summarized  by  Waste  Man¬ 
agement  Institute  (1991).  However,  most  of  these  programs  cover  only  near¬ 
shore  areas,  usually  extending  no  more  than  3  miles  offshore.  Of  the  ongoing 
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monitoring  efforts,  the  only  program  that  conies  close  to  fulfilling  the  require¬ 
ments  for  long-term  water  quality  monitoring  of  the  Bight  is  EPA's  New  Yorit 
Bight  Water  Quality  Monitoring  Program.  The  following  variables  are 
routinely  measured:  DO,  salinity,  temperature,  fecal  colifbrm  bacteria,  phyto¬ 
plankton  abundance  and  species  composition,  and  chlorophyll.  Additionally, 
sediment  and  benthic  samples  are  analyzed  for  viruses  and  other  pathogens  and 
heavy  metals.  Sampling  is  conducted  from  May  through  October.  The  EPA 
monitoring  program  contains  many  stations  along  foe  Long  Island  and  New 
Jersey  coasts  and  within  foe  Apex.  Several  transects  included  in  this  sampling 
program  extend  to  about  15  miles  offshore.  Additionally,  sediment  samples 
are  collected  in  a  transect  along  foe  Hudson  Canyon. 

If  long-term  monitoring  for  status  and  trends  of  hypoxia  and  eutrophication 
is  pursued,  then  the  following  recommendations  are  made.  These  recommen¬ 
dations  are  driven  by  foe  need  to  minimize  costs.  Thus,  a  small  number  of 
water  quality  variables  would  be  sampled  at  select  locations.  Any  new  moni¬ 
toring  program  should  be  conducted  in  conjunction  with  foe  existing  EPA  New 
York  Bight  Water  Quality  Monitoring  Program  to  reduce  costs.  Total  nitrogen 
and  phosphorus  should  be  added  to  foe  EPA  program  at  a  few  select  stations 
representative  of  the  New  York/New  Jersey  nearshore  and  Apex  regions.  Total 
N  and  P,  chlorophyll,  and  DO  should  also  be  monitored  at  new  stations  repre¬ 
sentative  of  the  midshelf  and  shelf  regions  (i.e.,  foe  six  regions  shown  in 
Plate  23).  Thus,  foe  sampling  program  would  basically  consist  of  three  tran¬ 
sects  (i.e..  New  Jersey  coast,  Apex/Hudson  Canyon,  and  Long  Island  coast) 
extending  from  nearshore  to  the  shelf  with  a  minimum  of  three  stations  along 
each  transect.  Each  station  should  be  sampled  at  near  surface  and  near  bottom 
for  total  N  and  P,  chlorophyll,  DO,  salinity,  and  temperature,  as  a  minimum. 
Additionally,  information  on  phytoplankton  abundance  and  species  composition 
is  desirable. 

Chlorophyll  has  been  monitored  through  remote  sensing.  This  approach 
requires  some  surface  measurements  for  calibration.  NOAA  is  one  Federal 
agency  that  has  used  remote  sensing  within  foe  Bight  for  chlorophyll  analyses. 
Remote  sensing  could  be  designed  into  a  future  monitoring  program  to  reduce 
costs  while  increasing  coverage. 

Before  a  new  long-term  water  quality  monitoring  program  is  developed, 
various  agencies  and  other  participating  parties  should  convene  to  scope  such 
an  effort,  to  avoid  duplication  while  trying  to  meet  foe  needs  of  all  partici¬ 
pants.  Several  other  related,  ongoing  efforts  could  influence  the  direction  of 
any  future  monitoring  program.  For  example,  NOAA  and  EPA  are  conducting 
ecosystem  monitoring,  which  includes  nutrient  analyses,  along  foe  northeast 
shelf  as  part  of  foe  Marine  Monitoring,  Analysis,  and  Prediction  Program. 
Additionally,  these  two  agencies  are  jointly  evaluating  the  existing  monitoring 
programs  nationwide  to  recommend  future  direction  in  coastal  monitoring 
under  the  National  Coastal  Monitoring  Program,  which  resulted  from  foe  1992 
Act  of  the  same  name.  New,  related  programs  might  be  under  way  by  univer¬ 
sities  in  foe  region  (such  as  foe  Long-Term  Ecosystem  Observation  Program 
recently  initiated  by  Rutgers  University).  The  specifics  of  any  additional. 
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future  long-term  monitoring  for  eutrophication  would  depend  on  the  opportuni¬ 
ties  for  cooperation  with  these  other  monitoring  efforts. 
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Appendix  A 

Water  Quality  Model  Kinetic 
Formulations  and  Coefficients 


Draft  text  from  Chapter  4  (The  Water  Quality  Model)  of  the  Chesapeake 
Bay  Water  Quality  Model  (Cerco  and  Cole,  in  preparation)1  is  extracted  and 
reproduced  in  this  appendix  to  relate  the  water  quality  model  (WQM)  kinetic 
coefficients  selected  and  to  clarify  the  New  York  Bight  (NYB)  model  adapta¬ 
tions.  The  major  NYB  kinetic  adaptations  that  included  sediment  oxygen 
demand  (SOD)  and  ammonia  ( NH4-N)  sediment-water  column  fluxes  and  bot¬ 
tom  layer  denitrification  are  also  presented.  A  tabulation  of  the  kinetic 
coefficient  values  used  in  the  NYB  water  quality  model  is  presented  at  the 
concluron  of  this  appendix. 


Sediment  Oxygen  Demand  (SOD) 

Sediment  oxygen  demand  was  specified  as  temporally  and  spatially  constant 
(SOD^).  The  SOD  was  modulated  through  temperature  and  DO  concentra¬ 
tion  of  overlying  water.  The  SOD  formulation  was 

SOD  =  SOD ^  exp  (KSO  (T  -  TRSO ))  DOKKHSO  ♦  DO)  (Al) 


where 

KSO  =  effect  of  temperature  on  SOD,  °C‘ 

TRSO  =  reference  temperature  for  SOD,  °C 

KHSO  =  half-saturation  concentration  of  DO  for  SOD,  g  Oj  m  3 


1  See  References  at  the  end  of  the  main  text 
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The  formulation  for  sediment  NH»-N  release  was  similar  to  that  for  SOD: 
NH4  release  -  NHA  ^  exp  (KSNH4  (T  -  TRSNH4) 

KHSNH4RKHSNH4  *  DO) 


(A2) 


where 

KSNH4  ss  effect  of  temperature  on  NH«-N  sediment  release,  °C‘ 

TRSNH4  =  reference  temperature  for  sediment  NH4-N  release,  °C 

KHSNH4  »  half-saturation  concentration  DO  for  sediment  NH4-N  release, 
8  Oi  m'3 


Bottom  Layer  Denitrification 

Nitrate  transfer  across  the  sediment-water  interface  is  proportional  to  the 
concentration  difference  between  sediment  interstitial  water  and  foe  water 
column. 


Flux  =  RKSED  (NOjMd  -  NOitMJ 


(A3) 


In  systems  such  as  the  NYB,  NOJwmm  »  NOimi  so  that 
Flux  =  -RKSED  NOs 


(A4) 


RKSED  =  sediment-water  mass  transfer  rate,  m  day'1 

The  following  text  was  extracted  from  Chapter  4  of  the  Chesapeake  Bay 
Water  Quality  Model  (Cerco  and  Cole,  in  preparation)  and  modified  to  reflect 
the  NYB  application. 


State  Variables 


Algae  are  grouped  into  two  size-fractionated  classes:  net  plankton  (>20  pm) 
and  nanoplankton  (<20  pm). 
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Organic  cartoon 

Two  organic  carbon  sate  variables  are  considered:  dissolved  (DOC)  ad 
particulate  (POC). 


>■»« - - — 

imrogon 

Nitrogen  is  divided  into  organic  and  mineral  fractions.  Organic  nitrogen  • 

state  variables  are  dissolved  organic  nitrogen  (DON)  and  particulate  organic 
nitrogen  (PON).  Two  mineral  nitrogen  forms  are  considered:  ammonia  (NHJ 
and  nitrate  (NO,).  The  nitrate  state  variable  represents  the  sum  of  nitrate  plus 
nitrite. 


Dtsaolired  oxygen 

Dissolved  oxygen  is  a  central  component  of  die  water  quality  model. 


Salinity 

Salinity  is  a  conservative  tracer  that  allows  for  verification  of  die  transport 
component  of  die  model  and  facilitates  examination  of  conservation  of  mass. 
Salinity  also  influences  the  saturation  concentration  of  DO. 


Temperature 

Temperature  is  a  primary  determinant  of  die  rate  of  biochemical  reactions. 


Conservation  Off  Mass  Equation 


The  Chesapeake  Bay  water  quality  model,  CE-QUAL-ICM,  is  an  integrated 
compartment  box  model.  The  box  structure  was  selected  to  allow  maximum 
flexibility  for  adaptation  of  die  model  to  alternate  hydrodynamic  models. 

Boxes  in  CE-QUAL-ICM  correspond  to  cells  in  X-Y-Z  space  on  die  CH3D 
grid.  CE-QUAL-ICM  solves,  for  each  grid  cell  and  for  each  state  variable,  the 
conservation  of  mass  equation: 


«  7-1  7-1  «*y 


(A5) 


where 

Vi  =  volume  of  i*  compartment  (m3) 


a.,,,  ami,!  A  »*«---  « — -  »<-  j-«  m — ■»-  *  -  .1 — 

a  rtmm  KjumKf  moam  fsmsc 
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C,  «  concentration  in  i*  compartment  (g  ms) 

Qj  *  volumetric  flow  across  flow  face  j  of  i*  compartment  (ms  sec'1) 

C*;  *  concentration  in  flow  across  flow  face  j  (g  m'*> 

Aj  =  area  of  flow  face  ;  (m2) 

Dj  a  diffusion  coefficient  at  flow  face  j  (m2  sec'1) 

n  =  number  of  flow  faces  attached  to  i*  compartment 

S,  a  external  loads  and  kinetic  sources  and  sinks  in  i*  compartment 
(g  sec1) 

t,  x  =  temporal  and  spatial  coordinates,  respectively 

This  appendix  focuses  on  the  portion  of  the  temporal  derivative  attributed  to 
internal  sources  and  sinks  (kinetics).  Within  the  model,  kinetics  are  computed 
using  a  temporal  dimension  of  days,  for  consistency  with  reported  rate  coeffi¬ 
cients.  Kinetic  sources  and  sinks  are  converted  within  CE-QUAL-ICM  to  a 
dimension  of  seconds  before  employment  in  the  mass-balance  equation. 


Algae 

Algae  {day  a  central  role  in  the  carbon  and  nitrogen  cycles  that  comprise 
the  model  ecosystem.  Equations  governing  the  two  algal  groups  are  largely 
the  same.  Differences  among  groups  are  expressed  through  .  te  magnitudes  of 
parameters  in  the  equations.  In  describing  the  parameters,  the  letter  "x"  is  used 
as  a  "wild  card.”  As  needed  in  this  text,  and  within  the  model  code,  the  wild 
card  is  replaced  with  a  letter  that  indicates  a  specific  algal  group.  Characters 
that  indicate  each  algal  group  are: 

d  =  net  plankton 

g  =  nanoplankton 

Sources  and  sinks  of  algae  are: 

Sources:  Growth  (production) 

Sinks:  Settling 

Basal  metabolism 
Predation 

The  last  two  sinks  of  algae  are  grouped  under  the  heading  mortality.  The 
governing  equation  for  algal  biomass  is 
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Bx  ■ 


ft  -  Mx  -  WS* 


where 

Bx  k  biomass,  expressed  as  carbon  (g  C  m  5) 
ft  »  production  (day1) 

Mx  «  mortality  (day*1) 

WSx  s  settling  velocity  (m  day'1) 
z  8  vertical  coordinate  (m) 


(A6) 


prooucnon 


Production  by  phytoplankton  is  determined  by  the  availability  of  nutrients, 
by  the  intensity  of  light,  and  by  die  ambient  temperature.  The  effects  of  each 
are  considered  to  be  multiplicative: 

ft  =  PMxmmm  (A' 


where 

PMx  =  production  under  optimal  conditions  (day'1) 
f[N)  a  effect  of  suboptimal  nutrient  concentration  (0  <  f  <  1) 
f(I)  =  effect  of  suboptimal  illumination  (0  <  f  <  1) 
fiT)  8  effect  of  suboptimal  temperature  (0  <  f  <  1) 


Nutrients 

Carbon,  nitrogen,  and  phosphorus  are  the  primary  nutrients  required  for 
algal  growth.  Carbon  is  usually  available  in  excess,  and  nitrogen  was  assumed 
the  limiting  nutrient  in  the  NYB.  Therefore,  carbon  and  phosphorus  were  not 
considered.  The  effect  of  nitrogen  on  growth  is  described  by  die  formulation 
commonly  referred  to  as  "Monod  kinetics"  (Monod  1949)  in  which  growth  is 
dependent  upon  nitrogen  availability  at  low  nitrogen  concentrations  but  inde¬ 
pendent  of  nitrogen  at  high  concentrations: 
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KHNx  =  half-saturation  constant  for  nitrogen  uptake  (g  N  m'5) 


Light 

Algal  production  increases  as  a  function  of  ligta  intensity  until  an  optimal 
intensity  is  reached.  Beyond  the  optimal  intensity,  production  declines  as 
intensity  increases.  Steele’s  equation  (DiToro,  O’Connor,  and  Thomann  1971) 
describes  this  phenomenon: 


ad 


jr_ 

is 


exp 


1  - 


(A9) 


where 

/  =  illumination  rate  (Langleys  day1) 

IS  =  optimal  illumination  (Langleys  day'1) 

Steele’s  equation  describes  the  instantaneous  light  limitation  at  a  point  in 
space.  The  model,  however,  computes  processes  integrated  over  discrete  time 
intervals  and  aggregated  spatially  into  model  segments.  Therefore,  Steele’s 
equation  must  be  integrated  over  an  appropriate  time  interval  and  averaged 
over  the  thickness  of  each  model  segment  The  integration  interval  selected  is 
1  day.  This  interval  does  not  preclude  computation  steps  less  than  a  day  but 
frees  the  model  from  accounting  for  illumination  in  "real  time."  The  price 
paid  is  that  diurnal  fluctuations  in  algae  are  not  computed.  Assuming  light 
intensity  declines  exponentially  with  depth,  the  integrated,  averaged  form  of 
Steele’s  equation  is: 

AD  *  f«P  (a*)  -  exp  (a/)]  (AlO) 


ab  *  -  exp  [-KESS  {ZD  *  H)]  (AlOa) 


at  *  -  exp  (-KESS  ZD)  (AlOb) 
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lo  *  daily  illumination  at  water  surface  (Langleys  day*1) 

FD  s  fractional  daylength  (0  <  FD  <  1) 

KESS  ■  total  light  attenuation  coefficient  (m4) 

ZD  =  distance  from  water  surface  to  top  of  layer  (m) 

Light  attenuation  in  the  water  column  is  composed  of  two  fractions:  a  back¬ 
ground  value  dependent  on  water  color  and  concentration  of  suspended  parti¬ 
cles,  and  extinction  due  to  light  absorption  by  ambient  chlorophyll: 

KESS  *  KE  *  KECHL  _J _  T  Bx  (All) 

CCHLx 


where 

KESS  -  total  attenuation  affecting  algal  growth  (m4) 
KE  =  background  light  attenuation  (m4) 


KECHL  =  light  attenuation  coefficient  for  chlorophyll  a  (m2  mg4) 

CCHLx  =  carbon-to-chlorophyll  ratio  of  algal  group  x  (g  C  mg4  chi) 

Optimal  illumination  for  photosynthesis  depends  on  algal  taxonomy,  dura¬ 
tion  of  exposure,  temperature,  nutritional  status,  and  previous  acclimation. 
Variations  in  optimal  illumination  are  largely  due  to  adaptations  by  algae 
intended  to  maximize  production  in  a  variable  environment.  Steele  (1962) 
noted  the  result  of  adaptations  is  that  optimal  illumination  is  a  consistent  frac¬ 
tion  («  SO  percent)  of  daily  illumination.  Kroner  and  Nixon  (1978)  reported 
an  analogous  finding  that  maximum  algal  production  occurs  at  a  constant 
depth  (-1  m)  in  the  water  column.  Their  approach  is  adopted  here  so  that 
optimal  illumination  is  expressed: 

ISx  =  IOAVG  exp  (-KESS  DOPTx)  (A12> 


where 

IOAVG  =  adjusted  surface  illumination  (Langleys  da/1) 

DOPTx  -  depth  of  maximum  algal  production  (m) 

A  minimum,  ISMIN,  is  specified  for  optimal  illumination  so  that  algae  do 
not  thrive  at  extremely  low  light  levels.  The  time  required  for  algae  to  adapt 
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to  changes  in  illumination  is  recognized  by  computing  ISx  based  on  a  time- 
weighted  average  of  daily  illumination: 

IOAVG  -  0.7  lo  *02  11  *  0.1  12  (A13> 


where 

11  =  daily  illumination  1  day  preceding  model  day  (Langleys  day'1) 

12  -  daily  illumination  2  days  preceding  model  day  (Langleys  day'1) 


Temperature 

Algal  production  increases  as  a  function  of  temperature  until  an  optimum 
temperature  or  temperature  range  is  reached.  Above  die  optimum,  production 
declines  until  a  temperature  lethal  to  the  organisms  is  attained.  Numerous 
functional  representations  of  temperature  effects  are  available.  Inspection  of 
growth  versus  temperature  curves  indicates  that  a  function  similar  to  a 
Gaussian  probability  curve  provides  a  good  fit  to  observations: 

exp  l-KTGxl  ( T  -  TMxf]  when  T  £  TMx 
F{T)  =  (A14> 

exp  [~KTGx2  {TMx  -  Tf)  when  T  >  TMx 


where 

TMx  =  optimal  temperature  for  algal  growth  (°Q 
KTGxl  =  effect  of  temperature  below  TMx  on  growth  (°C2) 
KTGx2  =  effect  of  temperature  above  TMx  on  growth  (°C3) 

Settling 

Separate  settling  velocities  are  specified  for  each  algal  group. 


Mortality 

Mortality  of  phytoplankton  is  the  sum  of  two  processes: 

Mx  -  BMx  *  PRx  *  BGx  <A15> 
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BMx  ■  basal  metabolism  (day1) 

PRx  »  predation  by  zooplankton  (day1) 


As  employed  here,  basal  metabolism  is  the  sum  of  all  internal  processes 
that  decrease  algal  biomass.  A  portion  of  the  metabolism  is  respiration  and 
may  be  viewed  as  a  reversal  of  production.  In  respiration,  carbon  and  nutri¬ 
ents  are  returned  to  die  environment  accompanied  by  die  consumption  of  DO. 

A  second  internal  sink  of  biomass  is  the  excretion  (exudation)  of  dissolved 
organic  carbon. 

Respiration  cannot  proceed  in  the  absence  of  DO.  Basal  metabolism  cannot 
decrease  in  proportion  to  oxygen  availability,  however,  or  algae  would 
approach  immortality  under  anoxic  conditions.  To  solve  this  dilemma,  basal 
metabolism  is  considered  to  be  independent  of  DO  but  the  distribution  of 
metabolism  between  respiration  and  excretion  is  DO  dependent  When  oxygen 
is  freely  available,  respiration  is  a  large  fraction  of  the  total.  When  oxygen  is 
restricted,  excretion  becomes  dominant  Formulation  of  this  process  is  detailed 
in  the  text  that  describes  algal  effects  on  carbon  and  DO. 

Basal  metabolism  is  commonly  considered  to  be  an  exponentially  increasing 
function  of  temperature: 

BMx  =  BMRx  exp  [KTBx  ( T  -  TRx )]  (A16) 


where 

BMRx  =  metabolic  rate  at  TRx  (day*1) 

KTBx  =  effect  of  temperature  on  metabolism  (“C*1) 
TRx  =  reference  temperature  for  metabolism  (°Q 


KffMMTOOn 

The  predation  formulation  is  identical  to  that  for  basal  metabolism.  The 
difference  in  predation  and  basal  metabolism  lies  in  the  distribution  of  the  end 
products  of  these  processes. 

PRx  »  BPRx  exp  [KTBx  (T  -  TRx)]  (AlT> 
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where  BPRx  it  the  predation  rate  at  TRx  (day1). 


Effect  of  algae  on  organic  cartoon 

During  production  and  respiration,  algae  primarily  take  up  and  produce 
carbon  dioxide,  an  inorganic  form  not  considered  in  the  model.  A  small  frac¬ 
tion  of  basal  metabolism  is  exuded  as  DOC,  however,  aid  in  the  model  this 
fraction  increases  as  DO  becomes  deficient.  Algae  also  produce  organic  car¬ 
bon  through  the  effects  of  predation.  Zooplankton  take  up  and  redistribute 
algal  carbon  through  grazing,  assimilation,  respiration,  and  excretion.  Since 
zooplankton  are  not  included  in  the  model,  routing  of  algal  carbon  through 
zooplankton  is  simulated  by  empirical  distribution  coefficients.  The  effects  of 
algae  on  organic  carbon  are  expressed  as 


DOC 


|  JfCDx  ♦  (1  - 


FCDx) 


KHRx 


KHRx  *  DO 


BMx  *  FCDP  PRx\  Bx 


"  POC  =  FCLP  PRx  Bx 
dt 


(A18) 


(A19) 


•  • 


where 

FCDx  -  fraction  of  basal  metabolism  exuded  as  DOC 

KHRx  =  half-saturation  concentration  for  algal  DOC  excretion  (g  DO 

m'3) 

FCDP  =  fraction  of  dissolved  organic  carbon  produced  by  predation 
FCLP  =  fraction  of  particulate  carbon  produced  by  predation 
The  sum  of  the  two  predation  fractions  must  equal  unity. 


Effect  of  algos  on  nitrogen  • 

Algae  take  up  NH<  and  N03  during  production  and  release  NH4  and  organic 
nitrogen  through  mortality.  NO3  taken  up  is  reduced  internally  to  NH4  before 
synthesis  into  biomass  occurs  (Parsons,  Takahashi,  and  Hargrave  1984).  Trace 

concentrations  of  NH4  inhibit  reduction  of  NO,  so  that,  in  the  presence  of  NH4  q 

and  NO,,  NH*  is  utilized  first  The  "preference"  of  algae  for  NH*  can  be 
expressed  empirically  (Thom arm  and  Fitzpatrick  1982): 
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(A20) 


PNx  -  NH4 


[iBWfeV 


N03 

Nk4)(KHftx  *  N03) 


KHNx 

(NH4  *  NOS)  ( KHNx 


♦  AK>5)] 


where 

PNx  -  algal  preference  for  NH4  uptake  (0  <  PNx  <  1) 

Algal  biomass  is  expressed  in  units  of  carbon.  Algal  uptake  and  release  of 
nitrogen  are  quantified  through  a  proportionality  constant  that  represents  the 
average  ratio  of  nitrogen  to  carbon  in  algal  biomass.  As  with  carbon,  routing 
of  algal  nitrogen  through  zooplankton  is  represented  by  distribution  coeffi¬ 
cients.  The  effects  of  algae  on  the  nitrogen  state  variables  are  expressed  as 

NH4  =  (BMx  FNlx  *  PRx  FNIP  -  PNx  px)  ANCx  Bx  (A21) 


d 

Tt 


N03  =  ( PNx  -  1)  pX  ANCx  Bx 


(A22) 


°  DON  •  {BMx  FNDx  *  PRx  FNDP )  ANCx  Bx  (A23) 

dt 


d 

It 


PON  *  0 BMx  FNLx  +  PRx  FNLP )  ANCx  Bx 


(A24) 


where 

ANCx  =  nitrogen-to-carbon  ratio  of  algae  (g  N  g'1  Q 
FNlx  =  fraction  of  inorganic  nitrogen  produced  by  metabolism 
FNDx  =  fraction  of  dissolved  organic  nitrogen  produced  by  metabolism 
FNLx  -  fraction  of  particulate  nitrogen  produced  by  metabolism 
FNIP  =  fraction  of  inorganic  nitrogen  produced  by  predation 
FNDP  =  fraction  of  dissolved  organic  nitrogen  produced  by  predation 
FNLP  =  fraction  of  particulate  nitrogen  produced  by  predation 
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_ 


The  sums  of  the  metabolism  fractions  aid  the  predation  fractions  must  each 
equal  unity. 


Effect  of  algae  on  (flaaolved  oxygon 

Algae  produce  oxygen  during  photosynthesis  and  consume  oxygen  through 
respiration.  The  quantity  produced  depends  cm  die  form  of  nitrogen  taken  up. 
Since  oxygen  is  released  in  the  reduction  of  NO,,  more  oxygen  is  produced, 
per  unit  of  caibon  fixed,  when  N03  is  the  algal  nitrogen  source  than  when  NH» 
is  the  source.  Equations  describing  algal  uptake  of  carbon  and  nitrogen  and 
production  of  DO  are  (Morel  1983): 


106  CO,  +  16  N//4*  +  HfO;  +  106  Hfi 
-*  protoplasm  *  106  O,  +  15  H* 

106  CO,  +  16  NO,*  +  HJ>0;  +  122  Hfi  +  17  H' 


(A25) 


(A26) 


-*  protoplasm  +  138  O, 


When  NH4  is  the  nitrogen  source,  1  mole  oxygen  is  produced  per  mole  carbon 
dioxide  fixed.  When  NO,  is  the  nitrogen  source,  1.3  moles  oxygen  are  pro¬ 
duced  per  mole  carbon  dioxide  fixed. 

The  equation  that  describes  the  effect  of  algae  on  DO  in  the  model  is 

.A  DO  =  (1.3  -  0.3  PNx)  Px  - _ — _ BMx\ AOCR  Bx  (A27) 

Tt  KHRx  +  DO 


where  AOCR  is  the  dissolved  oxygen-to-carbon  ratio  in  respiration  (2.67  g  DO  • 

g 1  C). 

The  quantity  (1.3  -  0.3  PNx)  is  the  photosynthesis  ratio  and  expresses  the 

molar  quantity  of  oxygen  produced  per  mole  carbon  fixed.  The  photosynthesis 

ratio  approaches  unity  as  the  algal  preference  for  NH*  approaches  unity.  _ 


Organic  Carbon 

Carbon  fixed  by  primary  producers  undergoes  innumerable  transformations 
in  the  water  column.  Particulate  organic  carbon  is  converted  to  dissolved 
organic  carbon,  DOC  and  POC  are  incorporated  into  hetero trophic  biomass, 
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and  organic  carbon  is  respired  to  inorganic  carbon.  A  detailed  representation 
of  the  carbon  cycle  is  not  required  to  achieve  the  objectives  of  this  model  A 
reduced  system  is  conceived  consisting  of  the  following  elements: 

Phytoplankton  exudation 
Predation  on  phytoplankton 
Dissolution  of  particulate  carbon 
Heterotrophic  respiration 
Denitrification 
Settling 

In  the  reduced  system,  POC  ami  DOC  are  produced  by  predation  on  phyto¬ 
plankton.  The  POC  undergoes  first-order  dissolution  to  DOC.  DOC  produced 
by  phytoplankton  exudation,  predation,  and  dissolution  is  respired  or 
denitrified  at  a  first-order  rate  to  inorganic  carbon. 


Dissolution  and  respiration  rates 

Dissolution  and  respiration  rates  depend  on  the  availability  of  carbonaceous 
substrate  and  on  heterotrophic  activity.  Heterotrophic  activity  and  biomass 
have  been  correlated  with  algal  activity  and  biomass  across  a  wide  range  of 
natural  systems  (Bird  and  Kalff  1984;  Cole,  Findlay,  aiid  Pace  1988).  Conse¬ 
quently,  algal  biomass  can  be  incorporated  into  dissolution  and  respiration  rate 
formulations  as  a  surrogate  for  heterotrophic  activity.  The  correlation  between 
algae  and  heterotrophs  occurs  because  algae  produce  labile  carbon  that  fuels 
heterotrophic  activity.  Dissolution  and  respiration  processes  do  not  require  the 
presence  of  algae,  however,  and  may  be  fueled  entirely  by  allochthonous  car¬ 
bon  inputs.  Representation  of  dissolution  and  respiration  in  the  model  allows 
specification  of  algal-dependent  and  algal-independent  rates: 

KDOC  =  KDC  *  KDCALG  £  Bx  (A28) 

where 

KDOC  =  respiration  rate  of  dissolved  organic  carbon  (day1) 

KDC  =  minimum  respiration  rate  (day1) 

KDCALG  =  constant  that  relates  respiration  to  algal  biomass  (m3  g  '  C 
day1) 

KLPOC  =  KDC  +  KDCALG  £  Bx  (A29) 
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where 


i 


KLPOC  -  dissolution  rate  of  particulate  organic  carton  (day1) 

KLC  -  minimum  dissolution  rate  (day'1) 

KLCALG  =  constant  that  relates  dissolution  to  algal  biomass  (m3  g'1  C 
day1) 

Temperature  has  a  multiplicative  effect  on  dissolution  and  respiration  as 
expressed  by  Equation  A16  with  appropriate  changes  in  notation. 


Denitrification  • 

As  oxygen  is  depleted  from  natural  systems,  oxidation  of  organic  matter  is 
effected  by  the  reduction  of  alternate  oxidants  (in  standard  terminology 
referred  to  as  "alternate  electron  acceptors”)-  The  sequence  in  which  alternate 
acceptors  are  employed  is  determined  by  the  thermodynamics  of  oxidation- 
reduction  reactions.  The  first  substance  reduced  in  the  absence  of  oxygen  is 
nitrate.  One  representation  of  the  denitrification  reaction  is 

4 N03  +4 H*  +  5  CH20  =  2N2  +  7  HjO  +  SC02  (A30) 


The  model  representation  of  the  denitrification  reaction  differs  from  the 
balanced  redox  equation  since  only  NO{  and  CHtO  (as  DOC)  are  model  state 
variables.  The  model  representation  incorporates  a  temperature -dependent 
reaction  rate  and  notes  that  significant  denitrification  occurs  only  when  nitrate 
is  freely  available  and  dissolved  oxygen  is  depleted: 


DENIT 


KHODOC  N03 

KHODOC  +  DO  KHNDN  +  N03 


AANOX  KDOC  (A31) 


where 


DENIT  =  denitrification  rate  of  dissolved  organic  carbon  (day1) 

AANOX  =  ratio  of  denitrification  to  oxic  carbon  respiration  rate 

KHODOC  =  half-saturation  concentration  of  DO  required  for  oxic  respir¬ 
ation  (g  DO  m'3) 

KHNDN  =  half-saturation  concentration  of  nitrate  required  for  denitri¬ 
fication  (g  N  m  3) 


Temperature  has  a  multiplicative  effect  on  denitrification  as  expressed  by 
Equation  A16  with  appropriate  changes  in  notation. 
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± 


Mm _ _ £ 


E/p— fallow  that  sum  all  organic  cubon  sources  sod  sinks  in  the 
model  ecosystem. 


Dissolved  organic  carbon 

s  «*  -£{[rcD*  * <>  -  "**  -nm^Bs]tUx  * FCDP  tRx)  “ 

♦  KLPOC  fiD  LPOC  +  KRPOC  AT)  RPOC 
'  Tmgrnxs  «^^doc-  den mm  doc 


(A32) 


Paniculate  organic  carbon 


4-  POC  =  £  FCLP  Pta  Bx 

at 


(A33) 


-  ATLPOC  /(r)  LPOC  -  WSL  Jf-  LPOC 


where  WSL  is  the  settling  velocity  of  labile  particles  (m  day1). 


Nitrogen  • 

Nitrogen  undergoes  innumerable  transformations  in  the  water  column.  A 
reduced  system  is  conceived  that  includes  the  following  processes: 

Algal  production  and  metabolism  * 

Predation 

Hydrolysis  of  particulate  organic  nitrogen 

Mineralization  of  dissolved  organic  nitrogen 

Settling 

Nitrification  • 

Denitrification 

Effects  of  nitrogen  of  algal  production,  metabolism,  and  predation  have 
already  been  detailed.  Descriptions  of  hydrolysis,  mineralization,  nitrification 
and  denitrification  follow. 
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Hydrolysis  and  mineralization 

For  purposes  of  this  model,  hydrolysis  is  defined  as  die  process  by  which 
particulate  organic  nitrogen  (PON)  is  converted  to  the  dissolved  organic  form. 
Mineralization  is  defined  as  the  process  by  which  DON  is  converted  to  NH*. 
Conversion  of  PON  to  NH»  proceeds  through  the  sequence  of  hydrolysis  and 
mineralization.  Direct  mineralization  of  PON  does  not  occur.  Formulations 
for  hydrolysis  and  mineralization  are  based  on  the  following  assumptions: 

a.  Rates  of  hydrolysis  and  mineralization  ate  proportional  to  available 
substrate. 

b.  Rates  of  hydrolysis  and  mineralization  are  proportional  to  algal  biomass. 

c.  Hydrolysis  and  mineralization  are  accelerated  when  inorganic  nitrogen  is 
insufficient  to  supply  algal  demand. 

Assumption  a  states  that  hydrolysis  and  mineralization  cannot  proceed  in 
the  absence  of  PON  or  DON.  The  assumption  is  a  restatement  of  first-order 
kinetics.  Assumption  b  recognizes  that  nitrogen  transformation  rates  are 
influenced  by  the  biomass  of  hetero trophic  organisms  that  mediate  the  trans¬ 
formations.  While  bacteria  and  zooplankton  are  not  quantified  in  the  model, 
observations  in  numerous  systems  (Bird  and  Kalff  1984;  Cole,  Findlay,  and 
Pace  1988)  indicate  their  biomass  is  proportional  to  algal  biomass.  Conse¬ 
quently,  algal  biomass  is  an  indicator  of  heterotrophic  biomass  and  relation  of 
nitrogen  transformations  to  algal  biomass  is  appropriate.  Assumption  c  is 
based  largely  on  analogy  to  phosphorus  mineralization  for  which  low  phos¬ 
phate  concentration  stimulates  production  of  an  enzyme  that  mineralizes 
organic  phosphorus  to  phosphate. 

Formulations  for  the  mineralization  and  hydrolysis  rates,  consistent 
with  the  three  assumptions,  are 


# 


•  • 


KDON  =  KDN 


KHN 


KHN  +  NH4  *  N03 


KDNALG  £  Bx  (A34) 


where 


KDON  =  mineralization  rate  of  dissolved  organic  nitrogen  (day1) 

KDN  =  minimum  mineralization  rate  (day1) 

KHN  =  mean  half-saturation  constant  for  algal  nitrogen  uptake  (g  N 
m3) 

KDNALG  =  constant  that  relates  mineralization  to  algal  biomass 

(m3  g'1  C  day1) 
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®  •( 


KLPON  -  KLN  ♦ 


KHN 


KLNALG  Y,  Bx  (A35) 


KHN  *  NH%  ♦  N03 


where 

KLPON  =  hydrolysis  rate  of  labile  particulate  nitrogen  (day1) 

KLN  =  minimum  hydrolysis  rate  (day'1) 

KLNALG  =  constant  that  relates  hydrolysis  to  algal  biomass  (m3  g1 
C  day  '1) 

Mineralization  and  hydrolysis  rates  are  a  function  of  temperature  as  expressed 
by  Equation  A16  with  appropriate  changes  in  notation. 


Nitrification 

Nitrification  is  a  process  mediated  by  specialized  groups  of  autotrophic  • 

bacteria  that  obtain  energy  through  the  oxidation  of  ammonium  to  nitrite  and 
oxidation  of  nitrite  to  nitrate.  A  simplified  expression  for  complete  nitrifica¬ 
tion  is 


NHa*  *  202  -*  NO{  +  H20  *  2H  * 


(A36) 


The  equation  indicates  that  2  moles  of  oxygen  are  required  to  nitrify  1  mole 
of  ammonium  into  nitrate.  The  simplified  equation  is  not  strictly  true,  how¬ 
ever.  Cell  synthesis  by  nitrifying  bacteria  is  accomplished  by  the  fixation  of 
carbon  dioxide  so  that  less  than  2  moles  of  oxygen  are  consumed  per  mole 
ammonium  utilized  (Wezemak  and  Gannon  1968). 


In  this  study,  nitrification  is  modeled  as  a  function  of  available  ammonium, 
dissolved  oxygen,  and  temperature: 


NT  =  — ££ _ NH*  fiXf)  NTM 

KHONT  *  DO  KHNNT  +  NHA 


(A37) 


where 

NT  =  nitrification  rate  (g  N  m'3  day'1) 

NTM  -  maximum  nitrification  rate  at  optimal  temperature  (g  N  m 3 

day1) 

KHONT  -  half-saturation  constant  of  DO  required  for  nitrification  (g  DO 

m'3) 
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KHNNT  *  half-saturation  constant  of  NH4  required  for  nitrification  (g  N 
of3) 

The  optimal  temperature  for  nitrification  may  be  less  than  peak  tempera¬ 
tures  that  occur  in  coastal  waters.  To  allow  for  a  decrease  in  nitrification  at 
superoptimal  temperature,  the  effect  of  temperature  on  nitrification  is  modeled 
in  the  Gaussian  form  of  Equation  A14  with  appropriate  changes  in  notation. 

Effect  of  nitrification  on  ammonium 

4-  NH4  •  -NT  (A38) 

at 


Effect  of  nitrification  on  nitrate 

4.  N03  *  NT 
dt 


(A39) 


Effect  of  nitrification  on  dissolved  oxygen 
4-  DO  -  -AONT  NT 

dt 


(A40) 


J 


where  AONT  is  the  mass  DO  consumed  per  mass  NH4-N  nitrified  (4.33  g  "V.: 
g1  N). 


Effect  of  denitrification  on  nitrate 


The  effect  of  denitrification  on  DOC  has  been  detailed.  Denitrification 
removes  nitrate  from  the  system  in  stoichiometric  proportion  to  DOC  removal 
as  determined  by  Equation  A31: 


d 

■3? 


-ANDC  DENIT  flT)  DOC 


(A41) 


where  ANDC  is  the  mass  N03-N  reduced  per  mass  DOC  oxidized  (0.933  g  N 
g1  C). 


Ammonia 

The  equation  is  written  by  summing  all  previously  described  sources  and 
rinks: 


» . . . ft. 
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4-  NH4  -  Y,  CAMx  FNIx  *  PRx  FNIP  -  Wi  Bx)  ANCx  Bx 

w  MU 


(A42) 


♦  KDONfiT )  DON  -  NT 


Dissolved  organic  nitrogen 
DON  «  (BAfa  FNDx  *  PRx  FNDP)  ANCx  Bx 

Hi 


-  KDON  AT)  DON  *  KLPON  AT)  LPON 

♦  KRPON  AT)  RPON 

Paniculate  organic  nitrogen 
.  PON  *  Y,  (BMx  FNLx  +  PRx  FNLF>  ANCx  Bx 

Hi 

-  KLPON  AT)  LPON  -  WSL  LPON 


(A43) 


(A44) 


Nitrate 

NOJ  -  Y  (PNx  ~  D  Px  ANCx  Bx 

Hi 

♦  NT  -  ANDC  DENIT  ATP)  DOC 


(A45) 


Dissolved  Oxygen 

Sources  and  sinks  of  DO  in  the  water  column  include: 

Sources:  Algal  photosynthesis 
Atmospheric  reaention 

Sinks:  Algal  respiration 

Heterotrophic  respiration 
Nitrification 
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AddddClOO 

Retention  occur  only  in  the  model  segments  that  form  the  sir- water  inter¬ 
face.  The  effect  of  re  aeration  is 

*  DO  -  Z!L  (J DOs  -  DO)  (A46) 

at  H 


where 

KR  s  re  aeration  coefficient  (m  day'1) 

DOs  -  dissolved  oxygen  saturation  concentration  (g  DO  ms) 

The  surface  renewal  concept,  attributed  to  Danckweits  by  O'Connor  and 
Dobbins  (19S8),  indicates 

KR  -  (DL  R)m  <A47> 


where 

DL  =  molecular  diffusivity  of  oxygen  in  water  (»1.7  x  1(7*  m2  day'1) 

R  =  surface  renewal  rate 

Specification  of  the  surface  renewal  rate  is  die  fundamental  problem  in 
reaeration  theory.  O’Connor  and  Dobbins  (1958)  state  that,  in  isotropic  tur¬ 
bulence,  surface  renewal  can  be  approximated  as  the  ratio  of  stream  velocity  to 
depth.  The  renewal  rate  is  also  influenced  by  wind,  however  (O’Connor 
1983).  Influences  on  reaeration  of  temperature  (ASCE  1961)  and  salinity 
(Wen  et  al.  1984),  most  likely  effected  through  changes  in  diffusivity,  have 
been  measured.  No  single  theory  that  unites  all  these  factors  into  a  formula¬ 
tion  of  reaeration  in  an  estuary  is  available.  The  surface  renewal  concept  is 
retained  in  this  study  with  the  renewal  rate  treated  as  a  calibration  parameter. 

Saturation  dissolved  oxygen  concentration  is  computed  (Genet,  Smith,  and 
Sonnen  1974): 

DOs  =  14.5532  -  0.38217  T  *  0.0054258  T1 


-  CL  (1.665  x  10*  -  5.866  x  10*  T  <A48> 


♦  9.7%  x  10*  T1 

where  CL  is  the  chloride  concentration  (=  salinity/1.80655). 
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Summary  of  DO  aourcM  and  sinks 

The  complete  kinetics  tot  DO  are 


DO  • 


£  [<u  - 


0.3  PNx)  Px  - 


KHRx 


P°  BMx 1 
fix  *  dO  J 


AOCRBx 


-  AONT  NT  -  _— ££  —  /J)  AOCR  KDOC  DOC  (M9> 

RtidBGC  rZ5<5  7 


(DOj  -  DO) 


Salinity 

No  internal  sources  or  sinks  of  salinity  exist 


Temperature 

A  conservation  of  internal  energy  equation  can  be  written  analogous  to  the 
conservation  of  mass  equation.  The  only  source  or  sink  of  internal  energy 
considered  is  exchange  with  the  atmosphere.  Although  solar  radiation  can 
penetrate  several  meters  into  the  water  column,  radiation-induced  increases  in 
internal  energy  are  here  assigned  entirely  to  the  surface  model  layer. 

For  practical  purposes,  the  internal-energy  equation  can  be  written  as  a 
conservation  of  temperature  equation.  Change  of  temperature  due  to  atmo¬ 
spheric  exchange  is  considered  proportional  to  die  temperature  difference 
between  the  water  surface  and  a  theoretical  equilibrium  temperature  (Edinger, 
Brady,  and  Geyer  1965): 

!^-p£v(ra-,,  <A50) 


where 

TE  -  equilibrium  temperature  (C°) 

KT  =  heat  exchange  coefficient  (watt  m'2  °C‘) 

CP  -  specific  heat  of  water  (4,200  watt  sec  kg'1  °C') 
p  =  density  of  water  (1,000  kg  m  3) 
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Kinetic  Coefficient  Vaiuee  Ueed  in  the  NYB 
Application 


The  coefficients  libeled  “Hardwired”  were  coded,  and  the  coefficients 
labeled  “Input”  were  specified  in  input  files. 


Hardwired: 

KE 

as 

0.15 

light  attenuation  coefficient,  m'1 

PMd 

» 

2.0 

maximum  net  plankton  production,  da/1 

PMg 

S 

2.3 

maximum  nanoplankton  production,  da/1 

BMRd 

s 

0.2 

metabolic  rate  net  plankton  at  reference  temperature. 

BMRg 

= 

0.23 

day1 

metabolic  rate  nanoplankton  at  reference 
temperature,  da/1 

predation  rate  net  plankton  at  reference  temperature. 

BPRd 

0.0 

BPRg 

— 

0.08 

day1 

predation  rate  nanoplankton  at  reference 

WSd 

s 

5.0 

temperature,  day1 

settling  rate  net  plankton,  m  day1 

WSg 

s 

0.1 

settling  rate  nanoplankton,  m  day1 

WSL 

= 

2.0 

settling  rate  labile  particulates,  m  day1 

RKSED 

55 

0.5 

bottom  layer  denitrification  rate,  m  day1 

Input: 

KSO 

a 

0.0? 

effect  of  temperature  on  SOD,  °C1 

KSNH4 

a 

0.07 

effect  of  temperature  on  NH«-N  sediment  release, 

O^vl 

KHSO 

s 

2.0 

half-saturation  concentration  of  DO  for  SOD,  g  02 
m3 

half-saturation  concentration  DO  for  sediment 

KHSNH4 

3“ 

2.0 

TRSO 

s 

20.0 

NH*-N  release,  g  02  m*J 
reference  temperature  for  SOD,  °C 

TRSNH4 

= 

20.0 

reference  temperature  for  sediment  NH«-N  release, 
°c 

TMd 

= 

20.0 

V 

optimal  temperature  for  net  plankton  growth,  °C 

KTGdl 

S 

0.004 

effect  of  temperature  on  net  plankton  growth  below 

KTGd2 

s 

0.006 

TMd 

effect  of  temperature  on  net  plankton  growth  above 

TMg 

= 

20.0 

TMd 

optimal  temperature  for  nanoplankton  growth,  °C 

KTGgl 

5= 

0.008 

effect  of  temperature  on  nanoplankton  growth  below 

KTGg2 

s 

0.010 

TMg 

effect  of  temperature  on  nanoplankton  growth  above 

CCHLd 

s 

80.0 

TMg 

C/Chl  ratio  for  net  plankton 

CCHLg 

ss 

80.0 

C/Chl  ratio  for  nanoplankton 

Appandx  A  Water  Quality  Modal  Knofc  Coafficteote 


KTOfl 

KTGg2 

m 

p  p 

i  i 

effect  of  temperature  on  nanopUnkton  growth  below 
TMg 

effect  of  tempentuie  on  nanopUnkton  growth  above 

CCHLd 

a 

80.0 

TMg 

C/Chl  ratio  for  net  plankton 

CCHLg 

= 

80.0 

C/Chl  ratio  for  nanoplankton 
algal  self-shading,  nr  (g  Chi)*1 

KECHL 

* 

17.0 

DOPTd 

= 

1.0 

depth  maximum  net  plankton  production,  m 

DOPTg 

= 

1.0 

depth  maximum  nanopUnkton  production,  m 
minimum  optimal  illumination  rate,  Langleys  day'1 

ISMIN 

= 

40.0 

IOWT 

= 

0.7 

weight  current  illumination 

I1WT 

= 

0.2 

weight  illumination  one  preceding  day 

I2WT 

= 

0.1 

weight  illumination  two  preceding  days 

KHNd 

= 

0.01 

half-saturation  concentration  of  N  for  net  plankton  N 

KHNg 

_ 

0.001 

uptake,  g  N*3 

half-saturation  concentration  of  N  for  nanopUnkton 

TRd 

s 

20.0 

N  uptake,  g  N*3 

reference  temperature  for  base  net  plankton 

TRg 

= 

20.0 

metabolism,  °C 

reference  temperature  for  base  nanoplankton 

KTBd 

= 

0.069 

metabolism,  °C 

effect  of  temperature  on  base  net  plankton 

KTBg 

= 

0.069 

metabolism,  °C1 

effect  of  temperature  on  base  nanoplankton 

FCDd 

= 

0.1 

metabolism,  °C1 

fraction  net  plankton  metabolism  excreted  as  DOC 

FCDg 

= 

0.1 

fraction  nanoplankton  metabolism  excreted  as  DOC 

KHRd 

0.5 

half-saturation  concentration  of  DO  for  net  plankton 

KHRg 

= 

0.5 

DOC  excretion,  g  Qj  m*3 

half-saturation  concentration  of  DO  for  nanoplank¬ 

FCDP 

s 

0.1 

ton  DOC  excretion,  g  02  m*3 
fraction  algal  predation  excreted  as  DOC 

KHODOC 

= 

0.5 

half-saturation  concentration  of  DO  for  heterotrophic 

FCLP 

— 

0.45 

respiration,  g  02  m*3 

fraction  algal  predation  excreted  as  POC 

KDC 

= 

0.01 

minimum  mineralization  rate  for  DOC,  day*1 

KDCALG 

a 

0.0 

constant  that  relates  respiration  to  algal  biomass,  m3 

KLC 

= 

0.035 

g*1  C  day*1 

minimum  hydrolysis  rate  for  POC,  day*1 

KLCALG 

a 

0.0 

constant  that  relates  dissolution  to  algal  biomass,  m3 

KTMNL 

s 

0.069 

g*1  C  day"1 

effect  of  temperature  on  mineralization,  °C*1 

TRMNL 

a 

20.0 

reference  temperature  for  mineralization,  °C 

KTHDR 

a 

0.069 

effect  of  temperature  on  hydrolysis,  °C1 

TRHDR 

a 

20.0 

reference  temperature  for  hydrolysis,  °C 

AANOX 

a 

0.5 

ratio  of  denitrification  to  oxic  carbon  respiration  rate 

KHNDN 

a 

0.1 

half-saturation  concentration  of  N03-N  for  denitrifi¬ 

TMNT 

a 

30.0 

cation,  g  N  m*3 

optimal  temperature  for  nitrification,  °C 
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KTMTl 

= 

0.09 

effect  of  temperature  on  nitrification  below  TMNT, 

OC1 

KTMT2 

S 

0.09 

effect  of  temperature  cm  nitrification  above  TMNT, 

’C1 

• 

KHONT 

= 

1.0 

half-saturation  concentration  of  DO  for  nitrification. 
gC^m'3 

KHNNT 

1.0 

half-saturation  concentration  of  NH4-N  for  nitrifica¬ 
tion,  g  N  m'3 

NTM 

as 

0.05 

maximum  nitrification  rate,  day'1 

• 

KHNd 

s 

0.01 

half-saturation  concentration  of  N  for  net  plankton  N 
uptake,  g  N  m'3 

KHNg 

0.001 

half-saturation  concentration  of  N  for  nanoplankton 

N  uptake,  g  N  m'3 

FNId 

0.25 

fraction  NH4-N  produced  by  net  plankton 
metabolism 

• 

FNIg 

0.25 

fraction  NH4-N  produced  by  nanoplankton 
metabolism 

ANCd 

= 

0.167 

net  plankton  N/C  ratio 

ANCg 

= 

0.167 

nanoplankton  N/C  ratio 

FN1P 

s 

0.10 

fraction  NH4-N  produced  by  algal  predation 

• 

KDN 

= 

0.04 

mineralization  rate  for  DON,  day'1 

KLN 

= 

0.10 

hydrolysis  rate  for  PON,  day'1 

FNDd 

0.75 

fraction  of  DON  produced  by  net  plankton 
metabolism 

FNDg 

0.75 

fraction  of  DON  produced  by  nanoplankton 
metabolism 

• 

FNLd 

= 

0.0 

fraction  of  PON  produced  by  net  plankton 
metabolism 

FNLg 

= 

0.0 

fraction  of  PON  produced  by  nanoplankton 
metabolism 

FNDP 

= 

0.0 

fraction  of  DON  produced  by  algal  predation 

• 

FNLP 

a 

0.45 

fraction  of  PON  produced  by  algal  predation 

» 
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Concentrations 


ApewriixB  Condition* 


field 

Surface 

Middle 

tst&am 

Twparitur*,  *C  I 

1 

fial 

lal 

2 

7.0 

7.0 

3 

dal 

dal 

4 

7.0 

7.0 

5 

lal 

lal 

6 

lad 

lal 

7 

dal 

lal 

8 

lal 

lafi 

9 

lal 

lal 

10 

lal 

lal 

11 

lal 

lad 

12 

lal 

lad 

13 

lal 

lal 

14 

10.0 

la 1 

15 

lal 

lal 

Salinity,  •/• 

1 

26.1 

26.3 

2 

30.0 

30.0 

3 

30.4 

Hal 

4 

30.0 

30.0 

5 

30.7 

Hal 

6 

Hal 

Hal 

7 

32.0 

Hal 

8 

Hal 

Hal 

9 

31.6 

Hal 

10 

Hal 

32.8 

11 

32.4 

Hal 

12 

Hal 

Hal 

13 

32.8 

34.0 

14 

32.4 

33.9 

15 

32.7 

33.6 

Netplankton,  mg  chi  a'1 

1  0.5 

2.0 

0.2 

2 

0.5 

2.0 

0.2 

3 

0.5 

2.0 

0.2 

4 

0.5 

2.0 

0.2 

5 

0.5 

2.0 

0.2 

6 

0.5 

2.0 

0.2 

7 

0.5 

2.0 

0.2 

8 

0.5 

2.0 

0.2 

9 

lal 

lal 

dal 

10 

lal 

lal 

dal 

11 

0.5 

2.0 

0.2 

12 

0.5 

2.0 

0.2 

13 

0.5 

2.0 

0.2 

14 

lal 

0.8 

dal 

15 

lal 

lal 

0.6 

Motes  Underlined 

values  are 

aeasured 

▼•lues. 

B2 


■anoplankton,  mg  eh: 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

Dissolved  Organic  C. 
1 


Particulate  Organic 
1 


Srld 


Burfaca 


Middle 


Aesnnia,  g  N  m3 

1 

0.007 

2 

0.007 

3 

0.007 

4 

0.007 

S 

0.014 

6 

0.007 

7 

0.007 

8 

0.007 

9 

0.006 

10 

0.004 

11 

0.004 

12 

0.004 

13 

0.006 

14 

9x991 

IS 

9.904 

Nitrate,  9  M  m3 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

Dissolved  Organic  Nitrogen,  g  N  m'3 


1 

0.014 

2 

0.014 

3 

0.014 

4 

0.014 

5 

0.014 

6 

0.014 

7 

0.014 

8 

0.014 

9 

0.014 

10 

0.014 

11 

0.014 

12 

0.014 

13 

0.014 

14 

0.014 

IS 

0.014 

9.149 

0.014 

9.011 

0.014 

9.018 

9.006 

9,001 

9.001 

0-001 

Q.001 

0,001 

0-003 

9*993 

Q.001 

Q-993 


0.014 

0.014 

0.014 

0.014 

9*919 

0.014 

0.014 

0.014 

9,925 

0.028 

0.028 

0.028 

9*991 

0.028 

9,1.025 


0.140 

0.007 

0.004 

0.007 

0.004 

0.007 

0.004 

■9.QQ3 

0.0S3 

0.039 

0.022 

ShSll 

0.109 

9,191 

9-125 


0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 
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firid 

fiUElASft  Middle 

fis££sa 

Particulate  Organic  Nitrogen,  g  N  a'5 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.014 

0.014 

9 

0.014 

0.014 

10 

0.014 

0.014 

11 

0.014 

0.014 

12 

0.014 

0.014 

13 

0.014 

0.014 

14 

0.014 

0.014 

15 

0.014 

0.014 

Dissolved  Oxygen, 

g  O,  nr3 

1 

1x1 

10.1 

2 

10.3 

10.0 

3 

10.3 

1*1 

4 

10.3 

10.0 

5 

10.3 

10.0 

6 

10.3 

1x1 

7 

1*2 

8.6 

8 

10.6 

1*2 

9 

10.6 

1x1 

10 

10.6 

8.9 

11 

10.6 

1x1 

12 

10.4 

8.6 

13 

11*4 

1x1 

14 

10*4 

1*1 

15 

10.4 

8.6 
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mini!  •  mum  *  -minn 


Omtmi  214  (Aug) 
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(Sap) 

Surf 

2S& 

1 

20.0 

17.2 

2 

21.8 

10.4 

3 

21.8 

11.0 

4 

21.5 

10.0 

5 

21*1 

8.8 

6 

20.6 

11.6 

7 

18.6 

13.0 

8 

21.2 

21.2 

275 
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0.10 

Mg 

3 

0.20 

2.00 

0.10 

0.20 

0.75 

o.os 

••g 

4 

0*22 

0.10 

2.00 

0.10 

0.20 

0.75 

0.05 

«g 

5 

1.12  0.05 

0.32 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

6 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

M9 

7 

1.08  3.89  0.12 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

wg 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Data: 

153 

( Jun) 

Surf  Mid  Bot 

Surf 

Hid 

Bat 

Surf 

Hid 

Bot 

sag 

1 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

2 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

3 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

4 

0.50 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

5 

0.68  0.02 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

6 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

7 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Data: 

183 

(Jul) 

Surf  Mid  Bot 

Surf 

Hid 

fiS£ 

Surf 

Mid 

fist 

sag 

1 

0.20 

2.00 

0.10 

0.20 

0.75 

o.os 

sag 

2 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

3 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

4 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

5 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

6 

0.50 

5.00 

0.10 

0.20 

0.75 

0.05 

sag 

7 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

sag 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

C6 
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Data:  214  (Aug) 


Surf 

Mid 

Bfli 

d)UC£ 

Mid 

M 

asui 

Mid 

ftafc 

1 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

2 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

3 

0.10 

1.00 

0.05 

0.10 

0.25 

0.01 

4 

0.10 

1.00 

0.05 

0.10 

0.25 

0.01 

S 

0.10 

1.00 

0.05 

0.10 

0.25 

0.01 

6 

0.10 

1.00 

0.05 

0.10 

0.25 

0.01 

7 

0.20 

2.00 

0.10 

0.20 

0.75 

0.05 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

M 

III 

<*«P> 

Surf 

Mid 

flat 

fiUEi 

Mid 

Bsi 

Surf 

Mid 

1G& 

l 

QaQI 

laid 

lalfi 

0.79 

2*11 

la 22 

2 

0.04 

Stall 

0.24 

0.16 

0.27 

flail 

3 

0.04 

QjlII 

0.24 

0.16 

0.27 

flail 

4 

0.03 

SLl21 

0.01 

0.16 

0.27 

flail 

S 

flafll 

flail 

0.01 

0.10 

0.25 

flaflfi 

6 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

7 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

275 

(Oct) 

Surf 

Mid 

B2t 

Surf 

Mid 

fisfl 

Surf 

Mid 

fist 

1 

0.05 

0.50 

0.25 

0.15 

0.25 

0.10 

2 

0.05 

0.50 

0.25 

0.15 

0.25 

0.10 

3 

0.03 

0.25 

0.15 

0.15 

0.25 

0.10 

4 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

5 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

6 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

7 

0.02 

0.10 

0.01 

0.10 

0.25 

0.00 

8 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Appendix  C  Boundary  Conditions 


C7 


BlMOlvtd  Organic 

Carbon, 

g  C  a4 

Dates  092  (Apr) 

Surf 

Ififc 

■eg  1 

2.00 

1.00 

Mg  2 

2.00 

1.00 

Mg  3 

2.00 

1.00 

Mg  4 

2.00 

1.00 

Mg  5 

2.00 

1.00 

Mg  6 

2.00 

1.00 

Mg  7 

2.00 

1.00 

Mg  8 

0.00 

0.00 

Oates  122  (Nay) 

Surf 

fis& 

seg  1 

2.00 

1.00 

aeg  2 

2.00 

1.00 

eeg  3 

2.00 

1.00 

aeg  4 

2.00 

1.00 

aeg  5 

2.00 

1.00 

aeg  6 

2.00 

1.00 

aeg  7 

2.00 

1.00 

aeg  8 

0.00 

0.00 

Oates  153  (Jun) 

Surf 

B fit 

aeg  1 

2.00 

1.00 

aeg  2 

2.00 

1.00 

aeg  3 

2.00 

1.00 

aeg  4 

2.00 

1.00 

aeg  5 

2.00 

1.00 

aeg  6 

2.00 

1.00 

aeg  7 

2.00 

1.00 

aeg  8 

0.00 

0.00 

Dates  183  (Jul) 

Surf 

fiat 

aeg  1 

2.00 

1.00 

sag  2 

2.00 

1.00 

aeg  3 

2.00 

1.00 

aeg  4 

2.00 

1.00 

aeg  5 

2.00 

1.00 

aeg  6 

2.00 

1.00 

aeg  7 

2.00 

1.00 

aeg  8 

0.00 

0.00 

Cl 
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Dttai  214  (lag) 


2 

3 

4 

5 

6 
7 

a 


tuxl 

2*22 

2.42 

2*42 

2*22 

2.00 

2.00 


2.00 

0.00 


245  (S«p) 


1 

2 

3 

4 

5 

6 

7 

8 


Surf 

2*22 

2*22 

2*22 

2*22 


2.00 

2.00 

2.00 

0.00 


275  (Oct) 


1 

2 

3 

4 

5 

6 

7 

8 


Surf 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

0.00 


1*22 

1*22 

1*22 

1*22 


1.00 

1.00 


1.00 

0.00 


2 at 
1*22 
1*22 
1*22 
1*22 
1.00 


1.00 

1.00 

0.00 


BS2 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 


Part icu lata  Organic  Carbon, 

g  C 

Data: 

092 

(Apr) 

Surf 

fiat 

••5 

1 

0.10 

0.10 

sag 

2 

0.10 

0.10 

■•9 

3 

0.10 

0.10 

4 

0.10 

0.10 

sag 

S 

0.10 

0.10 

Mfl 

6 

aoi 

sag 

7 

0.12 

••9 

8 

0.00 

0.00 

Oatet 

122 

(May) 

Surf 

flat 

aag 

1 

0.10 

0.10 

aag 

2 

0.10 

0.10 

sag 

3 

0.10 

0.10 

sag 

4 

0.10 

0.10 

aag 

5 

0.10 

0.10 

sag 

6 

0.10 

0.10 

sag 

7 

0.10 

0.10 

aag 

8 

0.00 

0.00 

Data: 

153 

(Jun) 

Surf 

Bot 

sag 

1 

0.10 

0.10 

sag 

2 

0.10 

0.10 

sag 

3 

0.10 

0.10 

sag 

4 

0.10 

0.10 

sag 

5 

0.10 

0.10 

sag 

6 

0.10 

0.10 

sag 

7 

0.10 

0.10 

sag 

8 

0.00 

0.00 

Data: 

183 

(Jul) 

Surf 

Bot 

sag 

1 

0.10 

0.10 

sag 

2 

0.10 

0.10 

sag 

3 

0.10 

0.10 

sag 

4 

0.10 

0.10 

sag 

5 

0.10 

0.10 

sag 

6 

0.10 

0.10 

sag 

7 

0.10 

0.10 

sag 

8 

0.00 

0.00 

laonl*  -  g  M  ar* 

Datai  092  (Apr) 

auri 

tot 

**g  1 

0.004 

0.004 

Mg  2 

0.004 

0.004 

Mg  3 

0.004 

0.004 

Mg  4 

0.004 

0.004 

Mg  5 

0.012 

9,910 

Mg  6 

0.005 

0-023 

Mg  7 

0.005 

a,  sas 

Mg  8 

0.000 

0.000 

Datai  122  (May) 

Surf 

M 

Mg  1 

0.004 

0.004 

sag  2 

0.004 

0.004 

sag  3 

0.004 

0.004 

sag  4 

0.005 

0.005 

sag  S 

0.006 

0.008 

sag  6 

0.004 

0.004 

sag  7 

0.005 

0.022 

sag  8 

0.000 

0.000 

Datai  153  (Jun) 

Surf 

fisfc 

sag  1 

0.019 

0.033 

sag  2 

ttJUJ 

Q.024 

sag  3 

0.017 

0-023 

sag  4 

0.003 

0.015 

sag  5 

0.004 

a,. 094 

sag  6 

0.004 

0.004 

sag  7 

0.004 

0.004 

sag  8 

0.000 

0.000 

Datai  Jul  (183) 

Surf 

fisfc 

sag  1 

0.004 

0.004 

sag  2 

0.004 

0.004 

sag  3 

0.004 

0.004 

sag  4 

0.004 

0.004 

sag  5 

0.004 

0.004 

sag  6 

0.004 

0.004 

sag  7 

0.004 

0.004 

sag  8 

0.000 

0.000 
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Data i  Aug  (214) 


Mg  1 

Surf 

0.004 

0.004 

Mg  2 

0.004 

0.004 

Mg  3 

0.004 

0.004 

Mg  4 

0.004 

0.004 

Mg  5 

0.004 

0.004 

Mg  6 

0.004 

0.004 

Mg  7 

0.004 

0.004 

Mg  8 

0.000 

0.000 

Data:  245 

Mg  1 

(Sap) 

Surf 

0.004 

isfc 

0.004 

Mg  2 

0.004 

0.004 

Mg  3 

0.004 

0.004 

Mg  4 

0.004 

0.004 

Mg  5 

0.008 

0.004 

Mg  6 

0.004 

0.004 

Mg  7 

0.004 

0.004 

Mg  8 

0.000 

0.000 

Data:  275 

Mg  1 

(Oct) 

Surf 

0.004 

ast 

0.004 

aag  2 

0.004 

0.004 

sag  3 

0.004 

0.004 

sag  4 

0.004 

0.004 

sag  5 

0.004 

0.004 

sag  6 

0.004 

0.004 

sag  7 

0.C04 

0.004 

sag  8 

0.000 

0.000 

Appandbc  C  Boundary  Condition* 


Nitrite  plus  Nitrate,  g  N  ■" 
Data:  092  (Apr) 


Surf 

fist 

aag  1 

0.014 

0.042 

aeg  2 

0.014 

0.042 

eeg  3 

0.042 

0.280 

aeg  4 

0.042 

0.280 

aeg  5 

S.&39 

0.280 

aeg  6 

0.003 

Q-126 

aeg  7 

0.003 

0.035 

aag  8 

0.000 

0.000 

Data:  122  (May) 


Surf 

Bot 

aeg  1 

0.014 

0.042 

aeg  2 

0.014 

0.042 

aeg  3 

0.014 

0.280 

aeg  4 

0.021 

Q.420 

aeg  5 

Q-.Q06 

aeg  6 

0.014 

0.042 

aeg  7 

0-002 

0.046 

aeg  8 

0.000 

0.000 

Date:  1S3  (Jun) 


Surf 

Bot 

aeg  1 

0-022 

0.032 

aeg  2 

0.010 

0.042 

sag  3 

0.006 

QU22A 

aeg  4 

0.064 

ft,.  16.4 

aeg  5 

SLJB51 

0-336 

aeg  6 

0.014 

0.042 

aeg  7 

0.014 

0.042 

aeg  8 

0.000 

0.000 

Date:  183  (Jul) 


Surf 

Bot 

aeg  1 

0.007 

0.021 

aeg  2 

0.007 

0.021 

aeg  3 

0.007 

0.140 

aeg  4 

0.007 

0.140 

aeg  5 

0.007 

0.140 

aeg  6 

0.007 

0.021 

aeg  7 

0.007 

0.021 

aeg  8 

0.000 

0.000 
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Datat  214 

(Aug) 

Surf 

Bfifc 

ug  1 

0.003 

0.009 

Mg  2 

0.002 

0.062 

MO  3 

0.003 

Q,U1 

Mg  4 

0.003 

Mg  5 

0.003 

0.224 

Mg  6 

0.001 

0.042 

sag  7 

0.001 

0.042 

sag  8 

0.000 

0.000 

Data:  245 

(Sap) 

Surf 

fis& 

sag  1 

0.007 

0.014 

sag  2 

0.009 

0.046 

sag  3 

0.005 

0.226 

sag  4 

0.014 

0.224 

sag  5 

Q...91S 

0.273 

sag  6 

0.004 

0.042 

sag  7 

0.004 

0.042 

sag  8 

0.000 

0.000 

Data:  275 

(Oct) 

Surf 

Bot 

sag  1 

0.007 

0.014 

sag  2 

0.007 

0.014 

sag  3 

0.007 

0.140 

sag  4 

0.007 

0.140 

sag  5 

0.007 

0.140 

sag  6 

0.004 

0.042 

sag  7 

0.004 

0.042 

sag  8 

0.000 

0.000 

Appendix  C  Boundary  Condition* 


Diaaolved  Organic  Nitrogen, 

i»«J  s 

no  MMurwtnti  on  boundary  or  in  interior 

Date:  092  (Apr) 

ansi 

fififc 

aeg  1  0.014 

0.014 

aeg  2  0.014 

0.014 

aeg  3  0.014 

0.014 

aeg  4  0.014 

0.014 

aeg  S  0.014 

0.014 

sag  6  0 . 014 

0.014 

aeg  7  0.014 

0.014 

aeg  8  0.000 

0.000 

Oates  122  (May) 

Surf 

aeg  1  0.014 

0.014 

aeg  2  0.014 

0.014 

aeg  3  0.014 

0.014 

aeg  4  0.014 

0.014 

aeg  S  0.014 

0.014 

aeg  6  0.014 

0.014 

aeg  7  0.014 

0.014 

aeg  8  0 . 000 

0.000 

Oates  153  (Jun) 

Surf 

Bot 

aeg  1  0.014 

0.014 

aeg  2  0.014 

0.014 

aeg  3  0.014 

0.014 

aeg  4  0.014 

0.014 

aeg  5  0.014 

0.014 

aeg  6  0.014 

0.014 

aeg  7  0.014 

0.014 

aeg  8  0 . 000 

0.000 

Dates  183  (Jul) 

Surf 

Bot 

aeg  1  0.014 

0.014 

aeg  2  0.014 

0.014 

aeg  3  0.014 

0.014 

aeg  4  0.014 

0.014 

aeg  5  0.014 

0.014 

aeg  6  0.014 

0.014 

aeg  7  0.014 

0.014 

aeg  8  0 . 000 

0.000 
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Data*  214  (Aug) 


Surf 

Bfifc 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

245  (Sap) 

Surf 

fist 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

275  (Oct) 


Surf 

Bot 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

Appendix  C  Boundary  Condition* 


Particulate  Organic  Nitrogen, 

g  N  a'3 

Data: 

092 

(Apr) 

Surf 

fist 

aag 

1 

0.014 

0.014 

aag 

2 

0.014 

0.014 

••g 

3 

0.014 

0.014 

*•9 

4 

0.014 

0.014 

••9 

S 

0.010 

0.014 

mmg 

6 

0.010 

0.0  Ai 

•mg 

7 

0.013 

0.011 

sag 

8 

0.000 

0.000 

Data: 

122 

(May) 

Sacf. 

Bot 

sag 

1 

0.014 

0.014 

sag 

2 

0.014 

0.014 

sag 

3 

0.014 

0.014 

sag 

4 

0.014 

0.014 

sag 

5 

0.014 

0.014 

sag 

6 

0.014 

0.014 

sag 

7 

0.014 

0.014 

sag 

8 

0.000 

0.000 

Data: 

153 

( Jun) 

Surf 

Bot 

s eg 

1 

0.014 

0.014 

sag 

2 

0.014 

0.014 

sag 

3 

0.014 

0.014 

aag 

4 

0.014 

0.014 

sag 

5 

0.014 

0.014 

aag 

6 

0.014 

0.014 

sag 

7 

0.014 

0.014 

sag 

8 

0.000 

0.000 

Date: 

183 

(Jul) 

Surf 

Bot 

sag 

1 

0.014 

0.014 

sag 

2 

0.014 

0.014 

sag 

3 

0.014 

0.014 

sag 

4 

0.014 

0.014 

aag 

5 

0.014 

0.014 

sag 

6 

0.014 

0.014 

sag 

7 

0.014 

0.014 

sag 

8 

0.000 

0.000 
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Detet  214  (Aug) 


■•g 
■eg 
■eg  8 


fiiud 

M 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

S 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

24S 

(Sep) 

&U£f 

fi at 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

275 

(Oct) 

Surf 

Bot 

1 

0.014 

0.014 

2 

0.014 

0.014 

3 

0.014 

0.014 

4 

0.014 

0.014 

5 

0.014 

0.014 

6 

0.014 

0.014 

7 

0.014 

0.014 

8 

0.000 

0.000 

Appendix  C  Boundary  Conditions 


Diaaolvad  Osygm, 

g0,«4 

Datai  092  (Apr) 

lad 

lsfc 

Mg  1 

8.57 

7.57 

aag  2 

8.57 

7.57 

Mg  3 

8.57 

7.57 

Mg  4 

8.57 

7.57 

Mg  S 

8.57 

7.57 

Mg  6 

9.43 

7.57 

Mg  7 

9.43 

7.57 

Mg  8 

12.86 

12.86 

Data:  122  (Nay) 

Surf 

M 

aag  1 

2x22 

7.43 

aag  2 

8.57 

7.43 

aag  3 

8.57 

7.43 

aag  4 
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